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This volume is the habihtation dissertation of the author written at the Fac- 
ulty of Mathematics, Computer Science and Econometrics of the University 
of Zielona Gora. 

The aim of this work is to present, in self-contained form, results concern- 
ing fundamental and the most important questions related to linear stochastic 
Volterra equations of convolution type. The paper is devoted to study the exis- 
tence and some kind of regularity of solutions to stochastic Volterra equations 
in Hilbert space and the space of tempered distributions, as well. 

In recent years the theory of Volterra equations, particularly fractional 
ones, has undergone a big development. This is an emerging area of research 
with interesting mathematical questions and various important applications. 
The increasing interest in these equations comes from their applications to 
problems from physics and engeenering, particularly from viscoelasticity, heat 
conduction in materials with memory or electrodynamics with memory. 

The paper is divided into four chapters. The first two of them have an 
introductory character and provide deterministic and stochastic tools needed 
to study existence of solutions to the equations considered and their regu- 
larity. Chapter [3] is devoted to stochastic Volterra equations in a separable 
Hilbert space. In chapter |4] stochastic linear evolution equations in the space 
of distributions are studied. 

The work is based on some earlier papers of the author, however part of 
results is not yet published. 

I wish to express my sincere gratitude to prof. Jerzy Zabczyk for intro- 
ducing me into the world of stochastic Volterra equations and inspiring math- 
ematical discussions. I am particularly grateful to prof. Carlos Lizama for 
fruitful joint research, exclusively through internet connection, which results 
is a big part of this monograph. 
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Introduction 



The main results in brief 

In the paper, two general problems concerning linear stochastic evolution 
equations of convolution type are studied: existence of strong solutions to such 
stochastic Volterra equations in a Hilbert space and regularity of solutions to 
two classes of stochastic Volterra equations in spaces of distributions. 

First, we consider Volterra equations in a separable Hilbert space H of the 
form 

X{t)=X{Q)+ [ a{t-T)AX{T)dT+ [ 1'{T)dW{T), t>0, (0.1) 
Jq Jo 

where X{Q) ^ H, a ^ -^ioc(^+) ^ scalar kernel function and A is a closed 
linear unbounded operator with the dense domain D{A) equipped with the 
graph norm. In (|0.ip . tf'(i), t > is an appropriate stochastic process and 
M^(i)i t > is a cylindrical Wiener process; both processes are defined on a 
stochastic basis {O^T, {J-t)t>o, P)- 

Equation (jO.ip arises, in the deterministic case, in a variety of applications 
as model problems, see e.g. |Pr93j . |Ba01| and references therein. Well-known 
techniques like localization, perturbation and coordinate transformation al- 
low to transfer results for such problems to integro-differential equations. In 
these applications, the operator A typically is a differential operator acting 
in spatial variables, like the Laplacian, the Stokes operator, or the elastic- 
ity operator. The kernel function a{t) should be thought as a kernel like 
a(t) = e-'^V-'^/r{l3);r] > 0, /3 G (0,2). The stochastic approach to integral 
equations has been recently used due to the fact that the level of accuracy for 
a given model not always seems to be significantly changed with increasing 
model complexity. 

Our main results concerning (jO.ip . rely essentially on techniques using a 
strongly continuous family of operators S{t),t > 0, defined on the space H 
and called the resolvent. Hence, in what follows, we assume that the deter- 
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ministic version of equation (jO.ip is well-posed, that is, admits a resolvent 
S{t),t> 0. 

The stochastic Volterra equations of the form (jO.ip have been treated by 
many authors, see e.g. |CD96j . |CD97j . |CDP97j . |CD00j or [RSOO] . [RSOT] . In 
the first three papers stochastic Volterra equations are studied in connection 
with viscoelasticity and heat conduction in materials with memory. The pa- 
per due to Clement and Da Prato |CD96| is particularly significant because 
the authors have extended the well-known semigroup approach, applied to 
stochastic differential equations, to a subclass of the equation (jO.ip . In the 
next papers, weak and mild solutions to the equation (j0.1[) have been studied 
and some results like regularity of solutions or large deviations of equations 
have been given. The resolvent approach to stochastic Volterra equations, 
introduced in |CD96j . enables us to obtain new results in an elegant way, 
analogously like in semigroup case. In resolvent case, new difficulties appear 
because the family S{t), f > 0, in general do not create a semigroup. 

Our main results concerning the equations (jO.ip in the space H are the 
existence theorems of strong solutions to some classes of such equations. We 
provide existence of strong solutions to (jO.ip under different conditions on the 
kernel function. In some cases, we arrive at stochastic versions of fractional 
Volterra equations with corresponding a-times resolvent families Sa{t), t>0. 
The key role in our proofs is played by convergence of resolvent or a-times re- 
solvent families corresponding to deterministic versions of Volterra equations. 
These convergence theorems are resolvent analogies of the well-known Hille- 
Yosida theorem in semigroup case. Our convergence results generalize theo- 
rems due to Clement and Nohel [CN79| obtained for contraction semigroups. 
Having such resolvent analogies of the Hille-Yosida theorem, we proved that 
the stochastic convolutions arising in Volterra equations (|0.ip are strong so- 
lutions to ([OTTjl . 

In the remaining part of the paper we study two classes of equations of 
convolution type: the equation 

X{t,e)^Xoi0)+ [ b{t-T)AX{T,9)dT + Wr{t,9), < > 0, 6i e M'', (0.2) 

and the integro-differential stochastic equation with infinite delay 

X{t,e)^ f b{t-s)[AX{s,e) + Wr{s,0)]ds, t > 0, BeT'^, (0.3) 

where T'^ is a d-dimensional torus. The kernel function b is integrable on R+ 
and the class of operators A contains the Laplace operator and its fractional 
powers. In both equations (|0.2p and (|0.3p . Wr denotes a spatially homoge- 
neous Wiener process, which takes values in the space of tempered distribu- 
tions S'(R'^), d>l. 

Equations (|0.2p and (|0.3p are generalizations of stochastic heat and wave 
equations studied by many authors. Particularly, regularity problems of these 
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equations have attracted many authors. For an exhaustive bibliography we 
refer to [Le01| . 

In the paper, we consider existence of the solutions to (|0.2p and (|0.3p in 
the space S"(]R'') and next we derive conditions under which the solutions to 
(|0.2p and (|0.3p take values in function spaces. 

In the case of equation (j0.2p , the results have been obtained by using the 
resolvent operators corresponding to Volterra equations. The regularity results 
have been expressed in terms of the spectral measure /x and the covariance 
kernel F of the Wiener process Wr- Moreover, we give necessary and sufficient 
conditions for the existence of a limit measure to the equation (|0.2p . 

In the case of the equation (|0.3p . we study a particular case of weak so- 
lutions under the basis of an explicit representation of the solution to (10. Sp . 
The regularity results have been expressed in terms of the Fourier coefficients 
of the space covariance F of the process Wr- 



A guided tour through the paper 

Chapter [T] has a preliminary character. Its goal is to introduce the reader to 
the theory of deterministic Volterra equations in Banach space and to provide 
facts used in the paper. Section fl.il gives notations used in the paper and sec- 
tion 11.21 gives basic definitions connected with resolvent operators. Sections 
11.31 and 1 1 .41 contain definitions and facts concerning kernel functions, particu- 
larly regular ones for parabolic Volterra equations. Some ideas are illustrated 
by examples. Section [1.51 provides new results due to Karczewska and Lizama 
|KL07c) ■ that is, the approximation theorems. Theorems 1 1 . 1 91 and 1 1 . 201 not 
yet published. These results are resolvent analogies of the Hille-Yosida theo- 
rem in semigroup case and play the same role like the Hille-Yosida theorem 
does. These results concerning convergence of resolvents for the deterministic 
version of the equation (jO.ip in Banach space play the key role for existence 
theorems for strong solutions and they are used in Chapter [31 

Chapter [2| contains concepts and results from the infinite dimensional 
stochastic analysis recalled from well-known monographs |CP78j . |DZ92| and 
|GT95j . Among others, we recall an infinite dimensional versions of the Fubini 
theorem and the Ito formula. Additionally, we recall a construction, published 
in [Ka98 , of stochastic integral with respect to cylindrical Wiener process. 
The construction bases on the Ichikawa idea for the stochastic integral with 
respect to classical infinite dimensional Wiener process and it is an alternative 
to the construction given in [DZ92 . 

Chapter [3| contains the main results for stochastic Volterra equations in 
Hilbert space. In Section [3Tl1 we introduce the definitions of solutions to the 
equation (jO.ip and formulate auxiliary results being a framework for the main 
theorems. In Section [3T2| we prove existence of strong solutions for two classes 
of equation ()0.ip . Basing on convergence of resolvents obtained in Chapter [H 
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we can formulate Lemma 13.151 and Theorem 13.161 giving sufficient conditions 
under which stochastic convolution corresponding to (jO.ip is strong solution 
to the equation (jO.ip . Section [231 deals with the so-called fractional Volterra 
equations. First, we prove using other tools than in Section Fl.Sl approximation 
results, that is. Theorems 13.201 and I3.21[ for a-times resolvents corresponding 
to the fractional equations. Next, we prove the existence of strong solutions 
to fractional equations. These results are formulated in Lemma l3 . 291 and The- 
orem 13.301 In Section 13.41 we give several examples illustrating the class of 
equations fulfilling conditions of theorems providing existence of strong solu- 
tions. 

In Chapter|4]we study regularity of two classes of stochastic Volterra equa- 
tions in the space of tempered distributions. Section [4. II has an introductory 
character. It contais notions and facts concerning generalized and classical 
homogeneous Gaussian random fields needed in the sequel. In Section [4.21 we 
introduce the stochastic integral in the space of distributions. Then we formu- 
late Theorem l4.4l which characterizes the stochastic convolution corresponding 
to (|0.2p . The main regularity results obtained in Section 14.21 are collected in 
Theorems 14.61 14.71 and 14.81 These theorems give sufficient conditions under 
which solutions to the equation (|0.2p are function- valued and even continuous 
with respect to the space variable. These conditions are given in terms of the 
covariance kernel F of the Wiener process Wr and the spectral measure of 
Wr^ as well. The results obtained in this section are illustrated by several 
examples. Section 14.31 is a natural continuation of the previous one. In this 
section we give necessary and sufficient conditions for the existence of a limit 
measure to the equation (|0.2p and then we describe all limit measures to (|0.2I) . 
The main results of this section, that is Lemmas 14.131 - 14.141 and Theorems 
14.151 - I4.16i are in a sense analogous to those formulated in !DZ96[ Chapter 
6], obtained for semigroup case. 

Section 113 is devoted to regularity of solutions to the equation (|0.3p . Here 
we study a particular case of weak solutions basing on an explicit representa- 
tion of the solution to (|0.3p . We find the expression for the solution in terms 
of the kernel h and next we reduce the questions of regularity of solutions 
to problems arising in harmonic analysis. Our main results of this section, 
that is Theorem 14.211 and Proposition 14.221 provide necessary and sufficient 
conditions under which solutions to (j0.3p are function- valued. These condi- 
tions are given in terms of the Fourier coefficients of the covariance F of the 
Wiener process Wr- Additionally, Section [44l contains some corollaries which 
are consequences of the main results. 
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Deterministic Volterra equations 



This chapter contains notations and concepts concerning Volterra equations 
used throughout the monograph and coUects some resuhs necessary to make 
the work self-contained. The notations are standard and follow the book by 
Pruss [Pr93) . 

Section [Ol has a preliminary character. In section FOl we recall the defi- 
nition of the resolvent family to the deterministic Volterra equation and the 
concept of well-posedness connected with the resolvent. Kernel functions, par- 
ticularly fc-regular ones, are recalled in sections 11.31 and [1.41 The above men- 
tioned definitions and results are described and commented in detail in Priiss' 
monograph [Pr93] and appropriate references therein. 

Section 11.51 originates from [KL07c] . Theorems 11.191 and 11.201 yet non- 
published, are deterministic approximation theorems which play a key role 
for existence of strong solutions to stochastic Volterra equations. These ap- 
proximation theorems are resolvent analogies to Hille-Yosida's theorem for 
semigroup and play the same role as that theorem. 

1.1 Notations and preliminaries 

Let _B be a complex Banach space with the norm | • |. We consider in B the 
Volterra equation of the form 



In dTII), a e iioc(R+ 

; R) is a non-zero scalar kernel; for abbreviation we will 
write a G Xj'Qj.(R+). A is a closed unbounded linear operator in B with a 
dense domain D{A) and / is a continuous B-valued function. In the sequel we 
assume that the domain D{A) is equipped with the graph norm | • |yi of A, i.e. 
\x\a ■= \x\ + \Ax\ for X e D{A). Then {D{A), \ ■ \a) is a Banach space because 
A is closed (see e.g. [ENOOj ) and it is continuously and densely embedded into 




(1.1) 



{B,\-\). 
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By (j{A) and g{A) we shall denote spectrum and resolvent set of the op- 
erator A, respectively. 

The equation (|l.ip includes a big class of equations and is an abstract 
version of several deterministic problems, see e.g. Pr93 . For example, if a{t) = 
1 and / is a function of C^-class, the equation (|l.ip is equivalent to the Cauchy 
problem 



Analogously, in the case a{t) — t and / of C'^-class, the equation (|l.ip is 
equivalent to 

u{t) = Au{t)+f{t) with initial conditions u(0) = /(O) and u(0) = /(O). 

Several other examples of problems which lead to Volterra equation (jl.ip can 
be found in [Pr93| Section 5]. 

In the first three chapters of the monograph we shall use the abbreviation 



for the convolution of two functions g and h. 

In the paper we write a(t) for the kernel function a. The notation a(t) will 
mean the function and not the value of the function a at t. Such notation will 
allow to distinguish the function a{t) and the article a. 

If a function v G L^^^IM.^; B) is of exponential growth, i.e. 

e^'^^\v{t)\dt < +00 for some w £ R, we can define the Laplace transform 
of the function v 



In the whole paper we shall denote by v the Laplace transform of the function 

V. 

In the whole paper the operator norm will be denoted by 1 1 • 1 1 . 
1.2 Resolvents and well-posedness 

The concept of the resolvent is very important for the theory of linear Volterra 
equations. The so-called resolvent approach to the Volterra equation (|l.ip has 
been introduced many years ago, probably by Friedman and Shinbrot ^FSB^; 
recently the approach has been presented in detail in the great monograph 
by Priiss [Pr93| . The resolvent approach is a generalization of the semigroup 
approach. 

By S{t), t > 0, we shall denote the family of resolvent operators corre- 
sponding to the Volterra equation (|l.ip and defined as follows. 



u{t) ^ Au{t) + f{t) with m(0)==/(0). 





ReA > UJ. 
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Definition 1.1. (see e.g. fFr93f ) 

A family (>S'(t))t>o of bounded linear operators in the space B is called 
resolvent for il.l]) if the following conditions are satisfied: 

1. S{t) is strongly continuous on R+ and S{0) — I; 

2. S{t) commutes with the operator A, that is, S{t){D{A)) C D{A) and 
AS(t)x = S{t)Ax for all x € D{A) and t > 0; 

3. the following resolvent equation holds 

S{t)x^x + [ a{t-T)AS{T)xdT (1.2) 
for all X e D{A), t > 0. 

We shall assume that the equation p.ip is well-posed in the sense that 
(jl.ip admits the resolvent S{t), t > 0. (Precise definition of well-posedness 
is given in IPr93| ). That dcfintion is a direct extension of well-posedness of 
Cauchy problems. The lack of well-posedness of (|1.1[) leads to distribution 
resolvents, see e.g. |DI84) . 

Proposition 1.2. fJPr93[ Proposition 1.1]) The equation \l.l\l is well-posed 
if and only if \1.1]) admits a resolvent S(t). If this is the case then, in addition, 
the range R{a * S(t)) C D{A) for allt>0 and 

S{t)x = x + Af a{t~T)S{T)xdT for all xeH,t>0. (1.3) 
Jo 

Comment Let us emphasize that the resolvent S{t), i > 0, is determined 
by the operator A and the function a{t). Moreover, as a consequence of the 
strong continuity of S{t) we have for any T > 

snp\\S{t)\\<+oo . (1.4) 

t<T 

Suppose S{t) is the resolvent for p.ip and let G '^{^) be an eigenvalue 
of A with eigenvector x 0. Then 

S{t)x = s{t;n)x, t>0, (1.5) 

where s{t; fj.) is the solution of the one-dimensional Volterra equation 

s{t;n)+fi[ a{t- T)s{T;^i)dT t>0. (1.6) 

Jo 

By WjQ^(R+;i?) we denote the Sobolev space of order (l,p) of Bochner 
locally p-integrable functions acting from into the space B, see e.g. [ENOOj . 

Definition 1.3. A resolvent S{t), for the equation iLl]) . is called diffe- 
rentiable if S{-)x G Wi^^{R+;B) for any x G DiA) and there 
exists a function ip G ijQj.(R+) such that \S{t)x\ < (p{t)\x\A a.e. on R+, 
for every x G D{A). 
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Similarly, if S{t) is differentiable then 

S{t)x ^ ljLr{t]ijL)x, i>0, (1.7) 
where r(i; /x) is the solution of the one-dimensional equation 

r{t;fj.)+fi[ a{t- T)r{T;n)dT = a{t), t>0. (1.8) 
Jo 

In some special cases the functions s{t; fi) and r(t] /i) may be found 
explicitely. For example, for a{t) = 1, we have s{t;fi) — r{t; ^) = e~^*', t > 
0, /i e C. For a{t) = t, we obtain s{t; fx) = cos{^^t), r{t] (x) = sm.{yfjit) / ,JJi^ 
t>0, /i e C. 

Definition 1.4. Suppose S{t), t > 0, is a resolvent for il.l]) . S{t) is called 
exponentially bounded if there are constants M > 1 and G M such that 

\\S{t)\\ < Me'^\ for all t > 0. 

(Af, w) is called a type of S{t). 

Let us note that in contrary to the case of semigroups, not every resolvent 
needs to be exponentially bounded even if the kernel function a{t) belongs to 
L^(K+). The Volterra equation version of the Hille-Yosida theorem (see e.g. 
|Pr93[ Theorem 1.3]) provides the class of equations that admit exponentially 
bounded resolvents. An important class of kernels providing such class of 
resolvents are a{t) = t°'~^ / r{a), a € (0, 2). For details, counterexamples and 
comments we refer to |DP93j . 

1.3 Kernel functions 

Two classes of kernel functions defined below play a prominent role in the 
theory of Volterra equations. 

Definition 1.5. A C°° -function a : (0, +oo) ^ M is called completely 
monotonic (-!)"«(") (t) > for allt>0, n<E N. 

Definition 1.6. We say that function a G L^{0,T) is completely positive 

on [0, r] if for any fx>0, the solutions of the convolution equations il.6]) and 
lil.8\) satisfy s{t] ^) > and r{t\^) > on [0,T], respectively. 

We recall that if a G Lj'Q(.(M+) is completely positive, then s{t;ii), the 
solution to (|1.6|) . is nonnegative and nonincreasing for any t > 0, /i > 0. In 
the consequence, one has < s{t; ji) < 1. This is a special case of the result 
due to LFr63 . 

Kernels with this property have been introduced by Clement and Nohel 
|CN79| . We note that the class of completely positive kernels appears quite 
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naturally in applications, particularly in the theory of viscoelasticity. Several 
properties and examples of such kernels appear in [Pr93[ Section 4.2]. 

^a — 1 

Examples 1. Let a{t) — — — — -, a > 0, where F is the gamma function. For 

r{a) 

a e (0, 1], function a{t) is completely monotonia and completely positive. 

2. Another example of completely positive function is a{t) — e~*, i > 0. 
An easy computation shows that then s{t; ^) = (1 + + Me^^^+^^*], for 

t,fi> 0. 

1.4 Parabolic equations and regular kernels 

This section is devoted to the so-called parabolic Volterra equations defined 
by Priiss [P?9T] . 

Let _B be a complex Banach space and 

^{uj,9) := {A e C : \arg{\-uj)\ < 6}. 

Definition 1.7. ffPfM Definition 2.1]) A resolvent S(t) for jfTl]) is called 
analytic, if the function S{-) : R-|- — > L{B) admits analytic extension to a 
sector ^{0,6q) for some < 9o < tt/2. An analytic resolvent S{t) is said 
to be o/ analyticity type (ciJo,^o) for each 6 < 9o and lu > luq there is 
M = M {uj, 6) such that 

||^WII<Me"^''^ z€Y.^O,0o). (1.9) 

Corollary 1.8. (JFEM Corollary 2.1]) Suppose S{t) is an analytic resolvent 
for 11. 1\) of analyticity type (wq, Oq). Then for each uj > loq and 9 < Oq there 
is M — AI(lu,9) such that 

||S'(")(<)|| < Mn!e'^*(i+")(at)-", t > 0, n G N, (1.10) 

where a — sin 9 . 

Analytic resolvents, the analog of analytic semigroups for Volterra equa- 
tions, have been introduced by Da Prato and lannelli 'DI80]. Analogously like 
in the theory of analytic semigroups, a characterization of analytic resolvents 
in terms of the spectrum of the operator A and the Laplace transform of the 
kernel function a{t) is possible. 

Theorem 1.9. (\Pr93[ Theorem 2.1]) Let A be a closed unbounded operator in 
B with dense domain D{A) and let a G ^[^^(R-i-) satisfy J* \a{t)\ e~^°-^dt < oo 
for some uja G K. Then \1.1\) admits an analytic resolvent S{t) of analyticity 
type {ll!q,9o) iff the following conditions hold: 

1) a{X) admits meromorphic extension to J2i^0i 9o + 7i'/2); 

2) a{X) ^ 0, and l/a(A) S g{A) for all X € J2i^o,9a + Tr/2); 
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3) For each > o-'rid 9 < 9q there is a constant C = C{uj,9) such that 
H(X) :— (l/a(A) — j4)^"'^/(A a(A)) satisfies estimate 

\\H{X)\\<C/\X-uj\ for all X e sum{Lj,9 + n/2). (1.11) 

Typical examples of the kernel functions a{t) and the operator A fulfilling 
conditions of Theorem 1 1.91 arc the following. 

Examples 1. Let kernels be a{t) = tl^-^/rijS), t > 0, where (3 £ (0, 2) and T 
denotes the gamma function. The pair {t'^"^ / r{(3), A) generates a bounded 
analytic resolvent iff g{A) D X](0,/37r/2) and WfJ-ifJ. ~ A)~^\\ < M for aU 

/xeE(0,/?^/2). 

2. An important class of kernels a{t) which satisfy the above conditions 
is the class of completely monotonic kernels. By [Pr93| Corollary 2.4] if addi- 
tionally a G C(0, oo) n L^(0, 1) and A generates an analytic semigroup T{t) 
such that \ \T{t)\ \ < M on S{0, 9) then (fTTTI) admits an analytic resolvent S{t) 
of type (0,0). 

Parabolic Volterra equations appear in a context of Volterra equations 
admiting analytical resolvents. 

Definition 1.10. Equation is caHerf parabolic, if the following condi- 

tions hold: 

1. a(A) 7^ 0, l/a(A) € q{A) for all ReA > 0. 

2. There is a constant M > 1 such that H{X) = (/ — a(X)A)^^ / X satisfies 
\\H{X)\\ < M/\X\ for all Re X > 0. 

From the resolvent point of view, the concept of parabolicity is between the 
bounded and the analytic resolvents: if (|l.ip admits an analytic resolvent S{t) 
then ()1.1|) is parabolic. On the other hand, if the equation (jl.ip is parabolic 
and the kernel function a{t) has some properties, like convexity, then the 
resolvent corresponding to p.ip has, roughly speaking, similar properties like 
analytic resolvent. 

Definition 1.11. Let a G Lj'Qj.(M+) be of subexponential growth and suppose 
a(A) ^ for all ReA > 0. The function a{t) is called sectorial with angle 
9 > (or merely (?-sectorialJ if \ arga(A)| < 9 for all ReA > 0. 

The standard situation leading to parabolic equations is provided by sectorial 
kernels and some closed linear densely defined operators A. 
The following criteria provide parabolic equations. 

Proposition 1.12. f \Pr93l Proposition 3.1]) Let a e Ll^^(R+) be 9-sectorial 
for some 9 < t:, suppose A is closed linear densely defined, such that g{A) D 
X;(0,6'), and\\{p~A)-^\\ < M/\n\ for all ^ G E(0'^)- 'Then fTT]) is parabolic. 
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The particular case is when A is the generator of a bounded analytic Co- 
semigroup and the function a{t) is 7r/2-sectorial. Because a{t) is 7r/2-sectorial 
if and only if a{t) is of positive type, we obtain the following class of parabolic 
equations. 

Corollary 1.13. ( \Pr93X Corollary 3.1]) Let a G L\^^{R+) he of subexponen- 
tial growth and of positive type, and let A generate a hounded analytical Cq - 
semigroup in B. Then il.l]) is parabolic. 

In the sequel we will need some regular kernels. 

Definition 1.14. Let a G ij'Q^(R+) be of subexponential growth and fc € N. 
The function a(t) is called fc-regular if there is a constant c > such that 
|A"a(")(A)| < c\a{X)\ for allReX > 0, < n < fc. 

Comment Any fc-regular kernel a{t), k > 1 has the property that a(A) has 
no zeros in the open right halfplane. 

We would like to emphasize that convolutions of fc-regular kernels are 
again fc-regular what follows from the product rule of convolutions. The inte- 
gration and differentiation preserve fc-regularity, as well. Unfortunately, sums 
and differences of fc-regular kernels need not be fc-regular. We may check it 
taking a{t) = 1 and b{t) = t^ . However, if a{t) and h{t) are fc-regular and 
I arga(A) — arg6(A)| < 6* < tt. Re A > then a{t) + h{t) is fc-regular. 

If a{t) is real-valued and 1-regular then a{t) is sectorial. The converse of 
this is not true. As the counterexample we can take a{t) = 1 for t G [0, 1], 
a{t) = for t > 1. 

Proposition 1.15. ( \Pr93\ Proposition 3.2]) Suppose a E Ll^^{M.-^-) is such 
that a(A) admits analytic extension to '^(0,y^), where ip > n/2, and there is 
9 G (0, oo) .such that \ arga(A)| < 9 for all A G ^(0, (p). Then a{t) is k-regular 
for every fc G N. 

So, nonnegative and nonicreasing kernels are in general not 1-regular but 
if the kernel is also convex, then it is 1-regular. 

Definition 1.16. Let a E Ll^^{M.-^-) and fc > 2. The function a{t) is called 
k -moiLotone if a E C''-'^{0,oo), a^"'^ (t) > for all t > 0, 0<n<k-2, 

and (— 1)'^"^ a'*^"^' (i) is nonincreasing and convex. 

By definition, a 2-monotone kernel a{t) is nonnegative, nonicreasing and con- 
vex, and a(t) is completely monotonic if and only if a{t) is fc-monotone for all 
fc > 2. 

Proposition 1.17. HPrQSl Proposition 3.3]) Suppose a E iJ^^(R+) is (fc + 1)- 
monotone, fc > 1. Then a{t) is k-regular and of positive type. 

Now, we recall the main theorem on resolvents for parabolic Volterra equa- 
tions. 
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Theorem 1.18. (\Fr93[ Theorem 3.1]) Let B he a Banach space, A a closed 
linear operator in B with dense domain D{A), a g Lj'Q^(K_|_). Assume U.lp 
is parabolic, and a{t) is k-regular, for some fc > 1. Then there is a resolvent 
S € C'^~^((0, oo); L(B)) for and there is a constant M > 1 such that 

estimates 

< A-f, for all t > 0, n<k-l, (1.12) 

and 

\\t'' S'^''-^\t)-s'' S^^-^\s)\\<M\t-s\[l+\og-^], 0<s<t<oo, (1.13) 

t — s 

are valid. 



1.5 Approximation theorems 

In this paper the following results contained in |KL07c| concerning conver- 
gence of resolvents for the equation p.ip in Banach space B will play the key 
role. They extend some results of Clement and Nohel obtained in |CN79| for 
contraction semigroups. Theorems 11.1 9lfTT20l and Proposition 11.211 are not yet 
published. 

Theorem 1.19. Let A be the generator of a Co-semigroup in B and suppose 
the kernel function a{t) is completely positive. Then (A, a) admits an exponen- 
tially bounded resolvent S{t). Moreover, there exist bounded operators An such 
that {An, a) admit resolvent families Sn{t) satisfying \ \Sn{t)\ \ < Me"""* {M > 
1, Wo > 0) for allt>0, n€ N, and 

Sn{t)x S{t)x as n ^ +oo (1-14) 

for all X B, t > 0. Additionally, the convergence is uniform in t on every 
compact subset o/R+. 

Proof. The first assertion follows directly from [Pr87[ Theorem 5] (see also 
|Pr93| Theorem 4.2]). Since A generates a Cg-semigroup T{t), t > 0, the 
resolvent set p{A) of A contains the ray [WjCxd) and 

\\RiKA)'^\\ < /'^ . forA>u;, fc G N, 
(A — ui)" 

where R{\, A) = {XI - A)-\ Xep{A). 
Define 

An := nAR{n, A) = n^R{n, A) - nl, n> w (1-15) 
the Yosida approximation of A. 
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Then 

oc 2k-i-k 
fc=0 

Hence, for n > 2w we obtain 



< Me""'*. (1.16) 



Taking into account the above estimate and the complete positivity of the 
kernel function a(i), we can follow the same steps as in |Pr87[ Theorem 5] to 
obtain that there exist constants M\ > and wi S K (independent of n, due 
to ([TTSD) such that 

||[i/n(A)]('^)|| < forA>w;i, 

A — Wi 

where H„(A) := (A — Aa(A)yl„)^^. Here and in the sequel the hat indicates 
the Laplace transform. Hence, the generation theorem for resolvent families 
implies that for each n e N, the pair {An, a) admits resolvent family 5„(i) 
such that 

||S'„(t)|| < Mie"'i* forallneN. (1.17) 
In particular, the Laplace transform ^^(A) exists and satisfies 

/>oo 

Sn{X)^Hn{X)= e-'^'Sn{t)dt, X>wi. 
Jo 

Now recall from semigroup theory that for all /i sufficiently large we have 

/•OO 



as well as, 

R{fi, A) = I e-^* T{t) dt . 



Since a(A) ^ as A ^ oo, we deduce that for all A sufficiently large, we 
have 

Aa(A) a(A) Aa(A) Jq 

and 

1 1 1 /-OO 
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Hence, from the identity 

and the fact that R{fi,An) R{fj.,A) as n — > oo for aU fj, sufficiently large 
(see e.g. |Pa83[ Lemma 7.3], we obtain that 

HniX)^H{X) asn^oo. (1.18) 

Finally, due to (|1.17p and (|1.18p we can use the Trotter-Kato theorem for 
resolvent families of operators (cf. [Li90l Theorem 2.1]) and the conclusion 
follows. □ 

Let us recall, e.g. from |Fa83| . that a family C{t), t > 0, of linear 
bounded operators on H is called cosine f Eunily if for every t, s > 0, t > s: 
C{t + s) + C[t - s) = 2C[t)C{s). 

Theorem 11.191 mav be reformulated in the following version. 

Theorem 1.20. Let A generate a cosine family T(t) in B such that\\T{t)\ \ < 
Me"* fort>0 holds, and suppose the kernel function a{t) is completely pos- 
itive. Then {A, a) admits an exponentially bounded resolvent S(t). Moreover, 
there exist bounded operators An such that (An , a) admit resolvent families 
Snit) satisfying ||5'„(i)|| < Me'""* (M > I, wq > 0) for allt>0, ne N, and 

Snit)x — > S{t)x as n ^ +oo 

for all X G B, t > 0. Additionally, the convergence is uniform in t on every 
compact subset o/M-i-. 

Remarks 

1. By [Pr93l Theorem 4.3] or [P?87l Theorem 6] Theorem [USD] holds also in 
two other cases: 

a) a{t) is a creep function with the function ai(t) log-convex; 

b) a = c-kc with some completely positive c £ Lj'Qj,(M+). 
(Let us recall the definition |Pr93|, Definition 4.4]: 

A function a : i-^- R is called a creep function if a{t) is nonnegative, 
nondecreasing, and concave. A creep function a{t) has a standard form 



a(t) = oq + Qoot + / ai{T)dT, t>0, 
Jo 



where oq = a(O-t-) > 0, aoo ~ Ihnt^oc a{t)/t — inff>oa(t)/t > 0, and 
ai{t) — a{t) — floo is nonnegative, nonincreasing, Imit^oc cii{t) — 0.) 

2. Other examples of the convergence (|1.14p for the resolvents are given, e.g., 
in |CN79| and [Fr69| . In the first paper, the operator A generates a linear con- 
tinuous contraction semigroup. In the second one, A belongs to some subclass 
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of sectorial operators and the kernel a(t) is an absolutely continuous function 
fulfilling some technical assumptions. 

Comment The above theorem gives a partial answer to the following open 
problem for a resolvent family S{t) generated by a pair {A, a): do exist 
bounded linear operators An generating resolvent families Sn{t) such that 
Sn{t)x S{t)x ? Note that in case a(t) = 1 the answer is yes, namely An are 
provided by the Hille-Yosida approximation of A and Sn (t) = e*'*" . 

The following result will be used in the sequel. 

Proposition 1.21. Let A, An and Sn{t) be given as in Theorem ] 1 . 1 9\ or The- 
orem \772Ul Then the operators Snit) commute with the operator A, for every 
n sufficiently large and t > 0. 

Proof. For each n sufficiently large the bounded operators An admit a re- 
solvent family S'„(t), so by the complex inversion formula for the Laplace 
transform we have 



where Fn is a simple closed rectifiable curve surrounding the spectrum of An 
in the positive sense. 

On the other hand, Hn{X) '■= (A — Xd{X)An) where A„ := nA{n — A)~^, 
so each A„ commutes with A on D{A) and then each -ff„(A) commutes with 
A, on D{A), too. 

Finally, because A is closed and all the following integrals are convergent 
(exist), for all n sufficiently large and x G D{A) we have 





□ 



2 



Probabilistic background 



In this chapter we recall from |CP78j . [Ic82) . |DZ92j and [GT95j basic and 
important concepts and results of the infinite dimensional stochastic analy- 
sis needed in the sequel. In particular, we present construction of stochastic 
integral with respect to a cylindrical Wiener process, published in |Ka98] . 

2.1 Notations and conventions 

Assume that {f2,J-',P) is a probability space equipped with an increasing 
family of cr-fields (Tt), t G I, where / = [0, T] or J = [0, +00), called filtration. 
We shall denote by Tt+ the intersection of all cr-fields !Fs, s > t. We say that 
filtration is normal if J-'q contains all sets B (z T with measure P{B) = and 
if Tt = Tt+ for any t € I, that is, the filtration is right continuous. 

In the paper we assume that filtration (J-t)t<^i is normal. This assump- 
tion enables to choose modifications of considered stochastic processes with 
required measurable properties. 

Let H and U be two separable Hilbert spaces. In the whole paper we 
write explicitely indexes indicating the appropriate space in norms | • |(.) and 
inner products (•,•)(.) • 

Definition 2.1. The H-valued process X(t), t E I , is adapted to the family 
{Ttjteij if for arbitrary t € I the random variable X{t) is Tt-measurable. 

Definition 2.2. The H-valued process X{t), t € [0, T], is progressively 
measurable if for every t € [0,T] the mapping [0,t] x — > iJ, {s,w) t-^ 
X^SjW) is B{[0,t]) X Tt- measurable. 

We will use the following well-known result, see e.g. [DZ92| . 

Proposition 2.3. (JDZM, Proposition 3.5]) LetX(t), t e [0, T], be a stochas- 
tically continuous and adapted process with values in H . Then X has a pro- 
gressively measurable modification. 
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By 7^00 we denote a cr-field of subsets of [0, +00) x SI defined as follows: Voo 
is the (T-field generated by sets of the form: (s, t] x F, where < s < t < +00, 
F £ J^s and {0} x F, when F £ J^q. The restriction of the cr-field Voo to 
[0,T]x n will be denoted by Vt- 

Definition 2.4. An arbitrary measurable mapping from ([0, +00) x fl,Voo) or 
([0, r] X Q,Vt) into {H; B{H)) is called a predictable process. 

Comment A predictable process is an adapted one. 

Proposition 2.5. (fDZM, Proposition 3.6]) Assume that X{t), t e [0,T], is 
an adapted and stochastically continuous process. Then the process X has a 
predictable version on [0, T] . 

By L{U,H), L{U) we denote spaces of linear bounded operators from U 
into H and in U , respectively. As previously, the operator norm is denoted by 
ll-ll- 

An important role will be played by the space of Hilbert-Schmidt opera- 
tors. Let us recall the following definition. 

Definition 2.6. (WKSTI or pZ9^ ) Assume that {cfc} C U and {fj} C H 
are orthonormal bases of U and H, respectively. A linear bounded operator 
T : U H is called Hilbert-Schmidt operator l^'^'cllf < -l-oo. 

Because 

k=l k=l j=l j = l 

where T* denotes the operator adjoint to T, then the definition of Hilbert- 

1 

Schmidt operator and the number HTH^/s = ^ do not depend 

on the basis {efc}, keN. Moreover \ \T\\hs = \ \T*\\hs- 

Additionally, L2{U,H) - the set of all Hilbert-Schmidt operators from U 
into H, endowed with the norm || • \ \hs defined above, is a separable Hilbert 
space. 



2.2 Classical infinite dimensional Wiener process 

Here we recall from jCP78j and |Ic82| the definition of Wiener process with 
values in a real separable Hilbert space U and the stochastic integral with 
respect to this process. 

Definition 2.7. Let Q : U U be a linear symmetric non-negative nuclear 
operator fTrQ < +00). A square integrable U -valued stochastic process W{t), 
t>0, defined on a probability space {f2,J^,{J^t)t>0yP), where Tt denote a- 
fields such that Tt d T for t < s, is called classical or genuine Wiener 

process with covariance operator Q if: 
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1. W{Q) = 0, 

2. Ew{t) = 0, GoY[W{t) - W{s)] ^(t- s)Q for all s,t>0, 

3. W has independent increments, that is W{s4) — W{s^) and W{s2) — W{si) 
are independent whenever < Si < S2 < S3 < S4, 

4-- W has continuous trajectories, 

5. W is adapted with respect to the filtration {J-t)t>o- 

If we choose to be the tr-field generated by {W^(s);0 < s < t}, then 
W{t) — W{s) is independent of JF, for a\lt>s from condition 3. of the above 
definition. Then E{W(t) - W{s)\T,} = E{W{t) ~ W{s)} = by condition 2. 
Hence, E{W(t)\J^s} = W{s) w.p.l and {W(t),Tt} is a martingale on [0, +00). 

We remark that an ahernative definition is to replace condition 4. of defini- 
tion l2.7l bv assuming that W{t) is Gaussian for all t > 0, see |CP78| for details. 

In the light of the above, Wiener process is Gaussian and has the following 
expansion (see e.g. |CP78[ Lemma 5.23]). Let {e^} C f/ be an orthonormal set 
of eigenvectors of Q with corresponding eigenvalues Q (so TtQ — X^i^i COj 
then 

00 

M^(i)-^A(Oe«, 

i=l 

where f3i are mutually independent real Wiener processes with E{0f{t)) = (it. 

Remark If W{t) is a Wiener process in U with covariance operator Q, then 
E\W{t) - W{s)\fp < {2n - iy.{t - s)"(TrQ)", where the equality holds for 
n = 1. 

The above type of structure of Wiener process will be used in the definition 
of the stochastic integral. 

For any Hilbert space H we denote by M{H) the space of all stochastic 
processes g : [0, T] x ^ L{U, H) such that 

E ^ \\9mliu,H)di^ < +^ 

and for all u E U, {g{t)u), t e [0, T] is an iJ-valued and JFf-adapted stochastic 
process. 

For each t e [0, T], the stochastic integral /J g{s)dW{s) £ H is defined for 
aU g G M{H) by 

/ g{s)dW{s) = lim V / g{s)em{s) 

Jo Jo 



in L^{[2;H). 

We shall show that the series in the above formula is convergent. 
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Let W^"'\t) = ^™ 1 ei(3i{t). Then, the integral 



is well defined for g G M (H) and additionally 



gis)dWis) 



in L^{fl]H). 



This convergence comes from the fact that the sequence 



Vm 



g{s)dW'^"'\s), mGN 



is Cauchy sequence in the space of square integrable random variables. Using 
properties of stochastic integrals with respect to Pi{s), for any m,n e N, 
m < n, we have: 

E {\yn - yrnlnj = (5(s)e»,g(s)ei)jj rfs (2.1) 

^ ( E c.u riiswii 



Hence, there exists a limit of the sequence (j/m) which defines the stochastic 
integral g{s)dW{s). 

The stochastic integral defined above has the following properties (see e.g. 
[k82]). 

Proposition 2.8. Let g £ M{H). Then 



g{t)dW{t) U ; 



2. E 



mdW{t) 



H 



E(Tr(.g(t)gg*(t)))dt 



E {TT{g*{t)g{t)Q)) dt < TrQ / E\git)\ldt 



3. E 



4. E 



sup 

te[o,T] 



sup 

tG[0,T] 



7(s)dW^(s) 



^{s)dW{s) 



2 " 








< 4E 


/ 9{s)dW{s) 






Jo 



2 J, 

<4:Ttq[ E\g{s)\jids 

H ^0 



H 



< 3E 



Tr(5(s)gff*(s))ds 
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(Since \ Jq g{s)dW{s)\'jf is a submartinagale, property 3 follows from Doob's 
inequality. Property 4 is also a consequence of a general inequality for mar- 
tingales.) 



2.3 Stochastic integral with respect to cyhndrical 
Wiener process 



The construction of the stochastic integral in section [2T2] required that Q was 
a nuclear operator. In some cases, this assumption seems to be artificial. For 
instance, all processes stationary with respect to space variable, have non- 
nuclear covariance operator. So, we shall extend the definition of the stochas- 
tic integral to the case of general bounded self-adjoint, non-negative operator 
Q on Hilbert space U. In this section we provide a construction, published in 
|Ka98| , of stochastic integral with respect to an infinite dimensional cylindri- 
cal Wiener process alternative to that given in [DZ92j . The construction is 
based on the stochastic integrals with respect to real-valued Wiener processes. 
The advantage of using of such a construction is that we can use basic results 
and arguments of the finite dimensional case. To avoid trivial complications 
we shall assume that Q is strictly positive, that is: Q is non-negative and 

QXy^OfOTXy^ 0. 

Let us introduce the subspace Uq of the space U defined by Uq — Q^{U) 
with the norm 



Q 



u e Uq. 

u 



Assume that [/i is an arbitrary Hilbert space such that U is continuously 
embedded into Ui and the embedding of Ua into C/i is a Hilbcrt-Schmidt 
operator. 

In particular, when Q — I, then Uq = U and the embedding of U into Ui is 
Hilbert-Schmidt operator. When Q is a nuclear operator, that is, Tr Q < -foo, 
then Uq = Q^{U) and we can take Ui = U. Because in this case Q ^ is Hilbert- 
Schmidt operator then the embedding Uq d U is Hilbert-Schmidt operator. 

We denote by L2 — L2{Uo,H) the space of Hilbert-Schmidt operators 
acting from Uq into H. 

Let us consider the norm of the operator ^6X2' 

00 00 

hM=l h,k=l 



where gj — y/Xjej, and {Aj}, {ej} are eigenvalues and eigenfunctions of the 
operator Q; 

{gj}, {ej} and {fj} are orthonormal bases of spaces Uq, U and H, respectively. 
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The space is a separable Hilbert space with the norm 
TiiiPQij*). 



Particular cases: 

1. If Q = / then Uq ~ U and the space becomes L2{U, H). 

2. When Q is a nuclear operator then L{U,H) C L2{Uo,H). Assume that 
K G L{U,H) and let us consider the operator = Klug, that is the 
restriction of operator K to the space Uq. Because Q is nuclear operator, 
then is Hilbert-Schmidt operator. So, the embedding J of the space Uq 
into U is Hilbert-Schmidt operator. We have to compute the norm 

of the operator : Uq^ H.We obtain HV'll^o = \\KJ\\lo = TtKJ{KJ)\ 
where J : Uq ^ U. 



Because J is Hilbert-Schmidt operator and K is linear bounded op- 
erator then, basing on the theory of Hilbert-Schmidt operators (e.g. 
|GV61j . Chapter I), KJ is Hilbert-Schmidt operator, too. Next, (KJ)* is 
Hilbert-Schmidt operator. In consequence, KJ{KJ)* is nuclear operator, 
so TiKJ{KJ)* < -l-oo. Hence, ip = K\iJo is Hilbert-Schmidt operator on 
the space Uq, that is K G L2{Uo, H). 

Although Propositions 12.91 and 12. lOl introduced below are known (see e.g. 
Proposition 4.11 in the monograph |DZ92] ). because of their importance we 
formulate them again. In both propositions, {.gj} denotes an orthonormal 
basis in Uq and {13 j} is a family of independent standard real- valued Wiener 
processes. 

Proposition 2.9. The formula 



t > 



(2.2) 



defines Wiener process in Ui with the covariance operator Qi such that 
TrQi < -l-oo. 



Proof. This comes from the fact that the series (|2.2I) is convergent in space 
L'^{n,J^,P;Ui). We have 



E 



E 



j—rn-\-l 



= E\ J2 9jPj{t), E 9Ait)\ =E {gMt),g,(3,{t)) 



Ui 



ij—m-{-l 



k—m-\'l 



Ui 



j=m+l 



= ^^| E (5,,5,)c/,/3|(i) =t E \9Ak^ 

ij—rrL-\-l j J— m-f-1 



n > m > 1. 
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From the assumption, the embedding J : Uq ^ Ui is Hilbert-Schmidt op- 
erator, then for the basis {gj}, complete and orthonormal in Uq, we have 
J2jLi l-^djlu^ < +00. Because Jg-j = gj for any g^ e Uq, then T.JLi \9j\u^ < 
-Kx) which means X)J=m+i \9j\ui ~^ ^ when m,n ^ 00. 

Conditions 1), 2), 3) and 5) of the definition of Wiener process are obvi- 
ously satisfied. The process defined by (|2.2p is Gaussian because /3j{t), j £ N, 
are independent Gaussian processes. By Kolmogorov test theorem (see e.g. 
|DZ92| . Theorem 3.3), trajectories of the process Wc{t) are continuous (con- 
dition 4) of the definition of Wiener process) because Wc{t) is Gaussian. 

Let Qi : Ui ^ Ui denote the covariance operator of the process Wc{t) 
defined by p.2|) . From the definition of covariance, for a,b G Ui we have: 



{Q^aM,^^E{a,W.{t)),^{h,WMK ^ E [Y,{a,g,),^{h,g,),^f3^^{t) 



= t^(a,5j)^^(6,5j)^^ y^gj{a,gj)^^,b 
j=i \i=i 



Hence Qio = t X^jli 5j (a,5j)[/^ • 

Because the covariance operator Qi is non-negative, then (by Proposition 
C.3 in |DZ92| ) Qi is a nuclear operator if and only if ^jOc/i < +00, 

where {/ij} is an orthonormal basis in Ui. 

From the above considerations 

00 00 00 

^{Qihj^hj)^^ < \gj\fj^ and then "^{Qikj , h^)^^ = TrQi < +00. 
j=i j=i j=i 

□ 

Proposition 2.10. For any a E U the process 

00 

{a,W,{t))^^Y.(^'9,)u(^j(t) (2-3) 
i=i 

is real-valued Wiener process and 

E (a, Wc{s))ij (6, Wc{t))^ ^(sAt) {Qa, b)^j for a,b£U. 

1 — i 

Additionally, ImQ^ = Uq and \u\iJo — Q^^u 

Proof. We shall prove that the series (|2.3p defining the process (a, Wc{t))jj is 
convergent in the space L^(J7;R). 
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Let us notice that the series (I2.3P is the sum of independent random vari- 
ables with zero mean. Then the series does converge in L^(i7;M) if and only 

if the following series J^'jLi ^ ((«> 9j)u l^i (^)) ^ converges. 
Because J is Hilbert-Schmidt operator, we obtain 



E 



2 

3lU 



Hence, the series 



does converge. Moreover, when t > s > 0, we have 



E ((a, WS))u (b, W^is))^) = E {{a, W,{t) - W,{s))^j (6, Wc{s))^j) 

+ E{{a,W,{s))^ib,WM)u) 
^E{{a,WrXs))^ib,W,{s))^) 



E 



.fe=i 



{b,9k)uPk{s) 



Let us introduce 



oc 



fe=i 



{b,gk)uf3k{t), for a,beU. 



Next, let S"^ and S*^ denote the partial sums of the series 5"^ and S'', respec- 
tively. From the above considerations the series 3°" and are convergent in 
L^{S1;R). Hence EiS^'S^) = limN ^oo E{S%S%). In fact, 



E\S%S% 



5a nb I TTil na nb na nb . nb na Qb na 



<E\S%\\S% 



b I q6 Qfl ( 

S''\ + E\S''\\S^ - S"] 



< {E\S%\^)^ {ElS^-S'l^y 

+ {Eis'fy {E\s% - s^fy 

because S"^ converges to S°' and 5^ converges to in quadratic mean. 
Additionally, E{SfjS'}f) — tJ^f^ii^'T gj)u(b, gj)u and when N +oo 



EiS'^S')=tJ2i<^,g,),ib,g,),. 



Let us notice that 



2.3 Stochastic integral with respect to cyhndrical Wiener process 

oo 

(Qia, b)^^ = E {a, W,{1%^ (b, W,{1))^^ = ^ (a, g,)^^ {b, g,)^^ 
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That gives Qi = J J*. In particular 

2 



(J*a,J*6)^^ = (JJ*a,6)j,^ 



Uo 



Qla 



Ui 



{Jra,a)jj^ = \ra\ 



Uo 



a £ Ui. 



(2.4) 



Having (|2.4|) . we can use theorems on images of linear operators, e.g. 
|DZ921 Appendix B.2, Proposition B.l (ii)]. 

By that proposition IvnQl — Im J. But for any j € N, and gj e Uq, 

Jgj — gj, that is Im J — Uq. Then IniQ^f" — Uq. 

_ 1 

Moreover, the operator G = ^ J is a bounded operator from Uq on Ui. 

_ 1 

From (|2.4p the adjoint operator G* — J*Qi is an isometry, so G is isometry, 



too. Thus 



□ 



Ui 



Ui 



In the case when Q is nuclear operator, is Hilbert-Schmidt operator. 
Taking Ui = U, the process Wc{t), i > 0, defined by (12. 2p is the classical 
Wiener process introduced in Definition 12.71 

Definition 2.11. The process Wc{t), t > 0, defined in i2.2\} . is called cylin- 
drical Wiener process on U whenTrQ = +oo. 

The stochastic integral with respect to cylindrical Wiener process is defined 
as follows. 

As we have already written above, the process Wc{t) defined by (|2.2p is 
a Wiener process in the space Ui with the covariance operator Qi such that 
TrQi < +00. Then the stochastic integral g{s)dWc{s) G H is well defined 
in Ui, where g{s) £ L{Ui,H). 

Let us notice that Ui is not uniquely determined. The space Ui can be 
an arbitrary Hilbert space such that U is continuously embedded into Ui and 
the embedding of Uq into C/i is a Hilbert-Schmidt operator. We would like to 
define the stochastic integral with respect to cylindrical Wiener proces Wc{t) 
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(given by (|2.2p ) in such a way that the integral is well defined on the space U 
and does not depend on the choice of the space Ui. 

We denote by A/'^(0, T; ij) the space of all stochastic processes 

[0,T]x L2{Uo,H) (2.5) 



such that 



\'P\\^r--=E[ I \\m\\Uuo,H)dt] (2-6) 



and for all u^Uq, ['P{t)u), t e [0, T], is an if -valued and -adapted stochas- 
tic process. 

The stochastic integral ^^<l>{s)dWc{s) € H with respect to cylindrical 
Wiener process, given by (|2.2p for any process € A/'^(0, T; i"): can be defined 
as the limit 



/ <?(s)(iW^e(s) = hm y / 
Jo ™-°°^Jo 



^>(s)ff,d/3,(s) (2.7) 



in h). 

Comment Before we prove that the stochastic integral given by the formula 
(|2.7p is well defined, let us recall properties of the operator Qi. From Propo- 
sition [2?9l cylindrical Wiener process Wc{t) given by (|2.2p has the covariance 
operator Qi : J7i — s- J7i, which is a nuclear operator in the space [/i, that is 

TrQi < +00. Next, basing on Proposition 12.101 Qi ■ Ui Uo, ImQl = Uq 



and lul 



Uo 



Qi 



for u e Uq- 

Moreover, from the above considerations and properties of the operator Qi 
we may deduce that L{Ui,H) C L2{Uo,H). This means that each operator 
# S L{Ui, H), that is linear and bounded from Ui into H, is Hilbert-Schmidt 
operator acting from Uq into H, that is € L2{Uo, H) when TrQi < +oo in 
Ui. This means that conditions and (EH) for the family N'^{Q,T; L!]) of 
integrands are natural assumptions for the stochastic integral given by (|2.7p . 

Now, we shall prove that the series from the right hand side of p.7p is 
convergent. 



Denote 



and 



Then, we have 



Z„r.^ f <I>{s)W^'^\s), te[0,T]. 
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2 

E \ I 'P(s)q,d0Js) for n > m > 1 



E \\Zn — Z„i\jj 



because from the assumption (|2.6 



^ / 1>{s)g,dP,{s) 
<E W{s)g,\lds 



H 



Then, the sequence {Z„i) is Cauchy sequence in the space of square- 
integrable random variables. So, the stochastic integral with respect to cylin- 
drical Wiener process given by (|2.7p is well defined. 

As we have already mentioned, the space Ui is not uniquely determined. 
Hence, the cylindrical Wiener proces Wdt) defined by (|2.2p is not uniquely 
determined either. 

Let us notice that the stochastic integral defined by (|2.7p does not depend 

on the choice of the space Ui. Firstly, in the formula (|2.7p there are not 

elements of the space Ui but only {gj}-basis of Uq. Additionally, in (|2.7p 

there are not eigenfunctions of the covariance operator Qi- Secondly, the class 

J\f'^{0,T; L2) of integrands does not depend on the choice of the space Ui 

- 

because (by Proposition l2.10p the spaces Ql{Ui) are identical for any spaces 

Ui: 

Qf : Ui^Ua and ImQf = C/q. 
2.3.1 Properties of the stochastic integral 

In this subsection we recall from |DZ92j some properties of stochastic integral 
used in the paper. 

Proposition 2.12. Assume that <P € A/'^(0, T; Lq). Then the stochastic inte- 
gral • W{t) /q ^s)dWis) is a continuous square integrable martingale 
and its quadratic variation is of the form « 'p9W{t) »— Q4,{s)ds, where 
Qs.(,s) = (<Z'(s)gV2)(^(,)Qi/2)*, s,te[0,T]. 

Proposition 2.13. If <P e Af'^{0,T;Ll), then E{<P»W{t)) = and 
E\<P»W{t)\% < +00, t e [0,T]. 

Proposition 2.14. Assume that <?i,^2 G -^^(0, T; Xq). Then the correlation 
operators V{s,t) := Cor(<?i • W{s),(I>2 • W{t)), s,t E [0,T] are given by the 
formula V{s,t) J^''\$2{r)Qi){<Pi{r)Qi)*dr. 
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Corollary 2.15. From the definition of the correlation operator we have 
E(<?i . Wis), $2 • W{t))H - E j;""* TT[i$2{r)Q^){-Pi{r)Qir]dr. 



2.4 The stochastic Fubini theorem and the Ito formula 

The below theorems are recalled directly from the book by Da Prato and 
Zabczyk |DZ92| . 

Assume that {Sl,T,{Tt)t>o,P) is a probability space, Ht ■= [0, T] x f2 
and recall that Vt is the cr- field defined in section [2TT1 that is Vt is the u-field 
generated by sets of the form: (s,t] x F, where 0<s<t<T,FE!Fs and 
{0} X F, when F e Tq. 

Let {E,£) be a measurable space and let 



<P : [t, Lu, x) — !■ uj, x) be a measurable mapping from 
{Qt E,Vt y-B{E)) into {Ll,B{Ll)), 



(2.8) 



where B{E) and ^(L") denote Borel cr- fields on E and L", respectively. Thus, 
in particular, for arbitrary x G E, <P{-, - ^x) is a predictable L2-valued process. 
Let in addition /i be a finite positive measure on (E,£). 

Theorem 2.16. (The stochastic Fubini theorem) 
Assume \2.8\) and that 



\\<P{-,-,x)\\t ^l{dx) < +00 



Then P-a. 



<^{t,x) dWit) 



^(dx) = 



(l>{t,x) fi{dx) 



dW{t) 



Assume that # is an L^-valned process stochastically integrable in [0, T], 
cj) is an iJ-valued predictable process Bochner integrable on [0,T], P-a.s., and 
X(0) is an .Fg-measurable _ff-valued random variable. Then the following pro- 
cess ^ ^ 

X{t)=X{0)+ [ (j){s)ds+ [ ^{s)dWis), t£[0,T], 
Jo Jo 

is well defined. Assume that a function F : [Q,T] x H ^ and its par- 
tial derivatives Ft, F^, Fxx, are uniformly continuous on bounded subsets of 
[0, T] X H. 
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Theorem 2.17. (The ltd formula) 
Under the above conditions 

F{t,X{t)) = F{0,X{0))+ I {F,{s,X{s)),<Pis)dW{s))H 

Jo 

{Ft{s,X{s)) + {F,{s,X{s)),^{s))H 
-^Tv[FUs,Xis)){<Pis)Q^){<P{s)Q-^r]}ds. 



holds P-a.s., for all t e [0,T]. 



3 



Stochastic Volterra equations in Hilbert space 



The aim of this chapter is to study some fundamental questions related to the 
linear convolution type stochastic Volterra equations of the form 



in a separable Hilbert space H . Particularly, we provide sufficient conditions 
for the existence of strong solutions to some classes of the equation (|3.1|1 . 
which is a stochastic version of the equation . 

Let (fl^T ^ {^t)t>Oj P) be a probability space. In (|3.ip . the kernel function 
a{t) and the operator A are the same as previously, X{0) is an i?- valued Tq- 
measurable random variable, is a cylindrical Wiener process on a separable 
Hilbert space U and is an appropriate process defined below. 

This chapter is organized as follows. In section [3TT] we give definitions of 
solutions to (|3.ip and some introductory results concerning stochastic convo- 
lution arising in (j3.1|) . Additionally, we show that under some conditions a 
weak solution to ()3.1|) is a mild solution and vice versa. These results have 
been recalled from Ka05j. 

Section 13.21 deals with strong solution to (|3.ip . We formulate sufficient 
conditions for a stochastic convolution to be a strong solution to (|3.ip . The 
above results come from the paper |KL07cj . not yet published. 

In section [331 based on [KL07d| . we study particular class of equations 
p.ip . that is, so-called fractional Volterra equations. We decided to consider 
that class of equations separately because of specific problems appearing dur- 
ing the study of such equations. First, we formulate the deterministic results 
which play the key role for stochastic results. We study in detail a-times resol- 
vent families corresponding to fractional Volterra equations. Next, we consider 
mild, weak and strong solutions to those equations. 

In the whole chapter we shall use the following Volterra Assumptions 




(3.1) 



(abbr. (VA)): 
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1. A : D{A) C H ^ H, is a. closed linear operator with the dense domain 
D(A) equipped with the graph norm | • 

2. a e Lj'Qj.(K+) is a scalar kernel; 

3. S{t), t > 0, are resolvent operators for the Volterra equation p.ip deter- 
mined by the operator A and the function a(t), t > 0. 

The domain D{A) is equipped with the graph norm defined as follows: 
■= (l^llf + l^'ill/)^ for h G D{A), where | • \h denotes a norm in H. 
Because H is a separable Hilbert space and ^ is a closed operator, the space 
{D{A), I • \d{A)) is a separable Hilbert space, too. 

t ^ 0; is a cylindrical Wiener process on U with the covariance 
operator Q and TrQ = -t-oo. 

By L2{Uq, H), as previously, we denote the set of all Hilbert-Schmidt 

operators acting from Uq into H, where Uq = Q^{U). 

For shortening, we introduce Probability Assumptions (abbr. (PA)): 

1. ^(0) is an i?- valued, ^o-measurable random variable; 

2. If" belongs to the space A^^(0, T; L2), where the finite interval [0, T] is fixed. 

3.1 Notions of solutions to stochastic Volterra equations 

In this section we introduce the definitions of solutions to the stochastic 
Volterra equation (|3.ip and then formulate some results, not yet published, 
setting a framework for further research. 

Definition 3.1. Assume that conditions (VA) and (PA) hold. An H -valued 
predictable process X{t), t €E [0,r], is said to he a strong solution to US. 

if X has a version such that P{X{t) G D{A)) ~ 1 /or almost all t £ [0, T]; for 
any t € [0, T] 

t 

|a(i - T)^X(r) |h dr < +00, P ~ a.s. (3.2) 
and for any t G [0, T] the equation I j*. 1\) holds P — a.s. 

Comment Because the integral /p ^{t) dW{T) is a continuous valued pro- 
cess then the above definition yields continuity of the strong solution. 

Let A* denotes the adjoint of the operator A, with dense domain D{A*) C 
H and the graph norm | • \d{A') defined as follows: |/i|d(a*) ■— (I^Ih + 
\A*h\jj)^, for h G D{A*). The space {D{A*), \ ■ \d{A')) is a separable Hilbert 
space. 
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Definition 3.2. Let conditions (VA) and (PA) hold. An H -valued predictable 
process X{t), t G [0,T], is said to be a weak solution to \3.1\] . if 
-P(/o \a{t-T)X{T)\HdT < +oo) = 1 and if for all £, e D(A*) and all t E [0,T] 
the following equation holds 

{X{t), Oh - (X(0), On + if a(t - t)X{t) dr, A* Oh 

+ {[ ^{T)dW{T),OH, P-a.s. 
Jo 

Definition 3.3. Assume that (VA) are satisfied andX{0) is an H -valued J-q- 
measurable random variable. An H -valued predictable process X{t), t £ [0, T], 
is said to be a mild solution to the .stochastic Volterra equation US. if 

e(^J^ \\S{t~T)1r{T)\\lodT^ < +00 for t<T (3.3) 

and, for arbitrary t G [0, T], 

X{t) = S{t)X{0)+ f S{t-T)<l^{T)dW{T), p-a.s. (3.4) 
Jo 

We will show that in some cases weak solution to the equation p.ip co- 
incides with mild solution to (|3.ip (see, subsection I3.1.2p . In consequence, 
having results for the convolution 

W'^{t):^f S{t-T)9{T)dW{T), te[0,T], (3.5) 
Jq 

where S{t) and are the same as in (|3.4p . we will obtain results for weak 
solution to (|3.ip . too. 

3.1.1 Introductory results 

In this section we collect some properties of the stochastic convolution of the 
form ^ 

W^{t):^ ( S{t^T)BdW{T) (3.6) 

"'0 

in the case when B S L{U, H) and is a cylindrical Wiener process. 

Lemma 3.4. Assume that the operators S{t), t > 0, and B are as above, 
S*{t),B* are their adjoints, and 

T pT 

\\S{r)B\\lodT^ Ti[S{t)BQB* S*{T)]dT <+oo. (3.7) 
^ Jo 

Then we have: 
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1. the process is Gaussian, mean-square continuous on [0,T] and then 
has a predictable version; 

2. 

Cov W^{t) =[ [SiT)BQB*S*{T)]dT, te[0,T]; (3.8) 

3. trajectories of the process are P-a.s. square integrahle on [0,T]. 

Proof. 1. Gaussianity of the process follows from the definition and 
properties of stochastic integral. Let us fix < t < < + < T. Then 

W^{t + h)-W^{t) ^ f [S{t + h-T)- S{t-T)]BdW{T) 

Jo 

j-t+h 

+ J S{t + h-T)BdW{T). 

Let us note that the above integrals are stochastically independent. Using 
the extension of the process W (mentioned in section [57T|) and properties 
of stochastic integral with respect to real Wiener processes (see e.g. 1^82) ). 
we have 

E\W''{t + h)-W''{t% = VAfc / \[S{t + h-T)-S{t~T)]Bek\jidT 

fc=l -^0 

+00 „t+h 

+ ^Afc/ \S{t + h - T)Bek\]i dT 

fc=i •'t 
:= Ii{t,h) + l2{t,h) . 

Then, invoking (|1.4p . the strong continuity of S{t) and the Lebesgue dom- 
inated convergence theorem, we can pass in Ii{t, h) with h under the 
sum and integral signs. Hence, we obtain Ii{t, h) ^ a.s h ^ 0. 
Observe that 

i-t+h 

h{tM=j^ \\S{t + h-T)BQmsdT , 

where || • \ \hs denotes the norm of Hilbert-Schmidt operator. By the con- 
dition p.7p we have 

/ \\S{t)BQ^\lsdt<+^, 
Ja 

what implies that lim/j^o hit, h) — 0. 

The proof for the case Q<t — h<t<T is similar. Existence of a 
predictable version is a consequence of the continuity and Proposition l2.5l 

2. The shape of the covariance p.8p follows from theory of stochastic integral, 
see e.g. |DZ92j . 
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3. From the definition (|3.6p and assumption p.7p we have the foUowing es- 
timate 

E 



/ \W^{r)\ldT= [ nW^{T)\ldT^ 
Jo J a 



T 



E 



S{t -r)BdW{r] 



dT^ / \\S{r)B\\lo drdT< +00. 

H Jo Jo ^ 



Hence, the function W^{-) may be regarded hke random variable with 
values in the space L^{Q, T; H). □ 

Now, we formulate an auxiliary result which will be used in the next sec- 
tion. 

Lemma 3.5. Let Volterra assumptions hold with the function a e 
VFi'i(R+). Assume that X is a weak solution to in the case when 

^{t) = B, where B S L(U,H) and trajectories of X are integrable P-a.s. on 
[0,T]. Then, for any function ^ e C^i[0,t]; D{A*)), t G [0,T], the following 
formula holds 

{X{t),at))H = {X{0),m)H+ f{{a^X){T)+a{Q)X{T),A*ar))HdT 

Jo 

+ / {^{T),BdW{T))H+ I {X{T),i{T))HdT, (3.9) 
Jo JO 

where dots above a and ^ mean time derivatives and * means the convolution. 

Proof. First, we consider functions of the form ^(t) := £^qip{t), t G [0, T], 
where ^o € D{A*) and (p G C^[0,T]. For simplicity we omit index h in the 
inner product. Let us denote F^g{t) := (X(t),^o), t G [0,T]. 
Using Ito's formula to the process F^„{t)ip{t), we have 

d[F^,{tMt)] = ipit)dF^„(t) + mFi„{t)dt, t e [o,r]. (3.10) 
Because X is weak solution to p.l|) . we have 



dF^o it) = {[ d{t - T)X{T)dT + a{0)X{t), A*(n)dt + {BdW{t), Co) 
Jo 

= ((a * X){t) + a{0)X{t),A*£,o)dt + {BdW{t), ^o)- (3.11) 
From (|XTU)) and (PUj) . we obtain 

F^„{tMt) = F^,iO)m + ! f{s){{a^X){s)+a{Q)X{s),A*^^)ds 

Jo 

+ I {^{s)BdW(s),^o)+ I ^{s){X(s),^o)ds 
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= {X{0),m)H+ f {id*X){s) + a{0)Xis),A*as))ds 
Jo 

+ f {BdW{s),as))+ [ {X{s),i{s))ds. 
Jo Jo 

Hence, we proved the formula (13. 9|) for functions ^ of the form ^(s) = ^o(p{s), 
s G [0, r]. Because such functions form a dense subspace in the space 
C^{[0,T];D{A*)), the proof is completed. □ 

3.1.2 Results in general case 

In this subsection we consider weak and mild solutions to the equation (j3.ip . 
First we study the stochastic convolution defined by (|3.5p . that is, 

W'^{t):= [ S{t-T)^{T)dW{T), ie[o,r]. 

^0 

Proposition 3.6. Assume that S{t), t > Q, are (as earlier) the resolvent 
operators corresponding to the Volterra equation Then, for arbitrary 

process ^ G 7V^(0,T;L2), the process t > Q, given by \3.5\l has a 

predictable version. 

Proof. Because proof of Proposition 13.61 is analogous to some schemes in the 
theory of stochastic integrals (see e.g. ^LS73, Chapter 4]) we provide only an 
outline of the proof. 

First, let us notice that the process S{t—T)'I'{T), where r S [0, t], belongs to 
Af^{0, T; L°), because if e Af^{0, T; L°). Then we may use the apparently well- 
known estimate (see e.g. Proposition 4.16 in |DZ92j ): for arbitrary a > 0, b > 
and t e [0, T] 

Pi\W''it)\H >a)<^ + P (^J^WSit - r)f (r)lliodr > b^ . (3.12) 

Because the resolvent operators S{t), t > 0, are uniformly bounded on com- 
pact itervals (see |Pr93| ). there exists a constant C > such that 1 15(^)1 1 < C 
for t G [0,T]. So, we have \\S{t - T)WiT)\\lo < C^\mT)\\lo, r G [0,T]. 
Then the estimate (I3.12p may be rewritten as 

Pi\W''it)\H>a)<^+p(^j\Mr)\\lodr>-^^ . (3.13) 

Let us consider predictability of the process in two steps. In the first 
step we assume that 'Z' is an elementary process understood in the sense given 
in section 4.2 in |DZ92] . In this case the process has a predictable version 
by Proposition 12. 51 
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In the second step ^ is an arbitrary process belonging to A/'^(0, T; Lj). 
Since elementary processes form a dense set in the space A/'^(0, T; Lj), there 
exists a sequence {^n ) of elementary processes such that for arbitrary c > 




(3.14) 



By the previous part of the proof the sequence W„ of convolutions 

W^{t) f S{t-T)Mr)dW{T) 
Jo 

converges in probability. Hence, it has a subsequence converging almost surely. 
This implies the predictability of the convolution W'^{t), t £ [0,T]. □ 



Proposition 3.7. Assume that W e 7V^(0, T;L^). Then the process W^it), 
t e [0,T], defined by iS. 5\) has square integrable trajectories. 



Proof. We have to prove that E {W^it^dt < 
and properties of stochastic integral 



From Fubini's theorem 



E 



S{t-T)^{T)dW{T) 



dt 



H 



\\S{t~T)W{T)\\lo drdt < M 



t 2 ■ 

S{t - T)>F{T)dW{T) 

H 



|<f(T)||2„ drdt < + 00 



dt 



(from boundness of operators S{t) and because '1^{t) are Hilbert-Schmidt). 



□ 



Proposition 3.8. Assume that a G BV{M.-^-), (VA) are satisfied and, addi- 
tionally S G C^(0, +00] L(H)). Let X be a predictable process with integrable 
trajectories. Assume that X has a version such that P{X{t) g D{A)) = 1 for 
almost all t e [0, T] and ( TO)) holds. If for any t G [0, T] and ^ e D{A*) 

{X{t),OH - (^(0), Oh + f\a{t - t)X{t), A*OHdT (3.15) 

Jo 

+ I {0^{T)dW{T))H, P^a.s., 
Jo 

then 

X{t) ^ S{t)X{0) + [ S(t-T)9{T)dW{T), te[0,T]. (3.16) 
Jo 
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Proof. For simplicity we omit index in the inner product. Since a G 
BV{M.+), we see, analogously like in Lemma [3.5[ that if (|3.15p is satisfied, 
then 

{Xit),at)) = (^(0),e(0))+ f\ia^ X){t) + aiO)XiT),A*aT))dT (3.17) 

Jo 

+ f (if(T)dW^(T),C(r)) + / (X(T),e(r))dT, P-a.s. 
Jo Jo 

holds for any $ e C\[0,t], D{A*)) for any t e [0,r]. 

Now, let us take C(t) := S'*(t-r)C with C e D{A*), t e [0, t]. The equation 
(|3.17p may be written like 

{X{t),S*{m = (X(0),5*(<)0+ l\{a*X){T)+a{0)X{T),A*S*{t-T)C}dT 



+ / {^{T)dW{T),S*{t~T)C)+ / (X(T),(5*(i-r)C)')dr, 



where derivative ()' in the last term is taken over t. 
Next, using S*{0) — /, we rewrite 



{X{t)X)^{S{t)X{Q),0+ {S{t-T)A 



d(T - a)X{a)da + a{0)X{T) 



+ / {S{t-T)rP{T)dW{T),0+ / (^(i-T)X(r),C)d- 



,OdT 

(3.18) 



To prove p.l6p it is enough to show that the sum of the first integral and the 
third one in the equation (|3.18|) gives zero. 

Because S € C^(0, +cx); L{H)) we can use properties of resolvent operators 
and the derivative S{t — r) with respect to r. Then 

/3 ( / S{t - T)X{T)dT, = (- / S{T)X{t - T)dT, 



a{T - s)AS{s)ds X{t~T)dT,C) 

d(r - s)AS{s)ds X{t - t^t- a{0)AS{T)X{t - T)d^ , 

= ((- [A{d * S) (r) * X] (t) - a{())A{S *X){t))X). 

Note that a G BV(M.^) and hence the convolution {a-kS){T) has sense (see 
[F?93l Section 1.6] or |AK89) V Since 



(a(0)AS'(t-T)X(r),Orfr = / (a(0)AS'(r)X(t - r), C)d- 
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and 



h 



{S{t~T)A 



a{T ~ (j)X {ij)da 



,OdT^ / (AS{t-T)id^X)iT)X)dT 



= (A(5*(a*X)(T))(i),C> = (^((^*a)(T)*X)W,C) for any C e D{A*), 

so /i = -/g. 

This means that p.l6p holds for any C G D{A*). Since is dense in 

H*, then (jXTC)) holds. □ 

Remark If (j3.ip is parabolic and the kernel a{t) is 3-monotone, understood 
in the sense defined by Priiss [Pr93| Section 3], then S G {0, +oo; L{H)), 
and a G BV{R+), respectively. 



Comment Proposition 13.81 shows that under particular conditions a weak 
solution to (|3.ip is a mild solution to the equation (|3.ip . 



Proposition 3.9. Let Volterra assumptions be satisfied. If 
'P E JV^{0,T; L2), then the stochastic convolution W'^ fulfills the equation 

mMi- 

Proof. Let us notice that the process has integrable trajectories. From 
the defintion of convolution ()3.5p . using Dirichlet's formula and stochastic 
Fubini's theorem, for any ^ e D{A*) we have 



{a{t-T)W'^{T),A*C)HdT= / {a{t-T) / S{t - cr)'F{<T)dW{a), A* ^ ndr 



a{t— t) S (t ~ a)dT 



^{<j)dW{a),A*£,)H 



t 

{A 



a{t — a—z)S{z)dz 



^{cT)dW{a),ClH ■ 



Next, using definition of convolution and the resolvent equation (|1.3p . as 
A(a -k S)(t — a)x = {S{t — a) — I)x, for x E H, wc can write 



{a{t-T)W^iT),A*OHdT^ { A[{a*S){t-a)]'F(a)dW{a),OH = 

Jq 

= {[ [Sit-a)-m'j)dWia),OH = 
Ja 



^{ S{t~a)^{a)dW{a),S)H -{ ^{a)dW{a), h ■ 

Jo JQ 

Hence, we obtained the following equation 



{W^it),OH= / {a{t^T)W^{T),A*OHdT+ / {^Mr)dWiT))H 
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for any ^ e I? (A*). □ 

Corollary 3.10. Let Volterra assumptions hold with a bounded operator 
A. If ^ belongs to M^{0,T; L°) then 



Comment The formula p.l9p says that the convolution W'^ is a strong 
solution to (|3.ip if the operator A is bounded. 

The below theorem is a consequence of the results obtained up to now. 

Theorem 3.11. Suppose that (VA) and (PA) hold. Then a strong solution 
(if exists) is always a weak solution of hS.l]) . If, additionally, assumptions of 
Provosition \3. 8\ are satisfied, a weak solution is a mild solution to the Volterra 
equation i3.1\) . Conversely, under conditions of Proposition\3^ a mild solu- 
tion X is also a weak solution to 113. 

Now, we provide two estimates for stochastic convolution p.Sp . 

Theorem 3.12. If e N'^{0,T]L%) then the following estimate holds 



where C is a constant and Mt = supf<y ||5'(t)||. 

Comment The estimate p.20p seems to be rather coarse. It comes directly 
from the definition of stochastic integral. Since (|3.20p reduces to the Davis 
inequality for martingales if S{-) = /, the constant C appears on the right 
hand side. Unfortunately, we can not use more refined tools, for instance Ito 
formula (see e.g. |Tu84| for Tubaro's estimate), because the process is not 
enough regular. 

The next result is a consequence of Theorem 13.121 
Theorem 3.13. Assume that <P e 7V2(0, T; L^). Then 



where a constant C{T) depends on T . 

Proof. From p.Sp and property of stochastic integral, writing out the Hilbert- 
Schmidt norm, we obain 




P-a.s. 



(3.19) 




(3.20) 



sup n\w''{t)\H)<c{T)mu-f^m-). 



t<T 
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n\W'^{t)\H)=Ei^\ S{t-T)1r{T)dW{T)\H 
<CE(^j\\Sit-TMT)\\lodT 



<CE[J^ \\S{t--r)\\'mr)\\iodTj < CMtE ^ mr)\\iodr 
< CMt (eJ^ ll'Z'Mlliodr^ ' < C{T) m\^.^o,T;LO) , 



where Mt is as above and C(T) = C Mr- □ 



3.2 Existence of strong solution 

In this section is a cyhndrical Wiener process, that is, TrQ — +oo and the 
spaces Uq, L2, A/'^(0,T;L2) ^re the same hke previously (see definitions in 
chapter The results from this section originate from IKL07cj and are not 
yet published. 

Let us recall the stochastic convolution introduced in (|3.5p 

iy*(t) := / S{t^T)^{T)dWiT), 



Jo 

where 'F belongs to the space A/'^(0, T; Lj)- In consequence, because resolvent 
operators S{t), t > 0, are bounded, then S{t - O^'l') ^ Af^iO, T; L^), too. 

In the sequel, by A<l>{t), t > 0, we will denote the composition of the 
operators <P{t) and A. 

We will use the following well-known result, where the operator A is, as 
previously, a closed linear operator with the dense domain D{A) equipped 
with the graph norm | • \d(A) and ^{t), t G [0,T] is an L2(?7o, i?)-predictable 
process. 

Proposition 3.14. (see e.g. \DZ9S\ Proposition 4-15]) 
If <Pit){Uo) C D{A), P - a.s. for all t e [0, T] and 



P 



then P 



£ wmwio dt<c^^^i, p (^£\\Am\\io dt<J^ = 1, 

^{t) dW{t) e D{A)^ = 1 and 

a[ <P{t)dW{t)= [ A<P{t)dW{t), P~a.s. 
Jo Jo 
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Let us recall assumptions of approximation theorems (Theorems 11.191 and 
ll.20p formulated for Hilbert space H: 

(ASl) The operator A is the generator of a Co-semigroup in H and the kernel 

function a[t) is completely positive. 
(AS2) A generates an exponentialy bounded cosine family in H and the 

function a{t) is completely positive (or fulfills one of two other cases listed 

in the remark 1 on pagell6|). 



Lemma 3.15. Let assumptions (VA) he satisfied. Suppose (ASl) or (AS2) 
hold. If^ and AW belong to M'^{0,T; L^) and in addition W{t){Uo) C D{A), 
P-a.s., then the following equality holds 

W'^{t)= [ a{t-T)AW'^{T)dT+ [ >F{T)dW{T), P-a.s. (3.21) 

Comment Let us emphasize that assumptions concerning the operators !f'(i), 
t > 0, particularly requirement that !Z'(i)(J7o) C D{A), P — a.s., are the same 
like in semigroup case, see e.g. |DZ921 Proposition 6.4]. 

Proof. Because formula p.l9p holds for any bounded operator, then it holds 
for the Yosida approximation An of the operator A, too, that is 

W^{t)^ I a{t-T)Ar,W^{T)dT+ f W{T)dWiT), 

Jo Jo 

where ^ 

W^{t) / Snit-TMr)dWir) 
Jo 

and ^ 

AnW^{t)=A^ I Sn{t-T)W{T)dW{T). 

Jo 

Recall that by assumption E J\f'^{0,T; L2). Because the operators S'„(t) 
are deterministic and bounded for any t E [0,T], n G N, then the operators 
Suit - ■)Wi-) belong to Af^{0, T; L^), too. In consequence, the difference 

<?„(t - •) := Sn{t - ■)W{-) - S(t - ■)>F{-) (3.22) 

belongs to 7V^(0, T; L^) for any t e [0, T] and n e N. This means that 




(3.23) 



for any t e [0,T]. 

Let us recall that the cylindrical Wiener process W{t), t > 0, can be 
written in the form 
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+00 



(3.24) 



where {gj} is an orthonormal basis of Uq and I3j{t) are independent real 
Wiener processes. From p.24p we have 



/ <l>n{t- t) dW{T) -^nit-r) g, dp, (r) 

Jo Jo 



(3.25) 



From (|3.23p . we obtain 



E 



< +00 



(3.26) 



for any t £ [0,r]. Next, from (|3.25p . properties of stochastic integral and 
((X^ we obtain for any t e [0, T], 



E 



<?„(i-r) dWir) 



E 



H 



E 



+00 „t 



< E 



+00 „f 

/ ^n{t^T)g,dP,{T 

+00 

E/ |<?„(T-r)5,|l,dT 



< 



< +00. 



By Theorem 11.191 or 11.201 the convergence (|1.14p of resolvent families is 
uniform in t on every compact subset of ]R-(-, particularly on the interval [0, T]. 
Now, we use (jl.l4p in the Hilbert space H, so (jl.l4p holds for every x E H. 
Then, for any fixed j, 



r I [5„ (T - r) - S{T ~ r)] 1r{r) g,\Ur — 
Jo 



(3.27) 



for n — > +00. Summing up our considerations, particularly using (|3.26p and 
(|3.27p we can write 



sup 1 
te[o,T] 



$nit~T)dW{T) 



= sup E 

H te[o,T] 



[Sn{t~T)~S{t~T)]^iT)dW{T) 



H 



< E 



+00 „T 

Y \[S„iT-T)-S{T-rmT)g,\Ur 



as n ^ +00 . 



Hence, by the Lebesgue dominated convergence theorem 
lim sup E\W'^{t)-W'^{t)\'^„^0. 



(3.28) 
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By assumption, If (t)(t/o) C ^(A), P~a.s. Because S{t){D{A)) C D{A), 
then S{t-T)^{T){Uo) C D{A), P - a.s., for any t e [0,t], t > 0. Hence, by 
Proposition E31 P(H/*(t) e D{A)) = 1. 

For any n G N, t > 0, we have 



where 



|A„W,f (i) - AW'^{t)\H < Nn,l{t) + NnMt), 



iV„4(i) := \Ar,W^{t) - A^W''it)\H, 

NnMt) ■■= \AnW^{t) - AW^{t)\H = |(A„ - A)W''{t)\H 



Then 



\A,,W^{t) - AW'^{t)\j, < Nl,{t) + 2N,,s{t)NnAt) + Nl,2{t) 

<ml.i{t)+Nl^{t)]. (3.29) 

Let us study the term Nn,i{t). Note that the unbounded operator A gen- 
erates a semigroup. Then we have for the Yosida approximation the foUowing 
properties: 



AnX = JnAx for any x £ D{A), sup || < oo 



(3.30) 



where AnX = nAR{n,A)x = AJnX for any x H, with J„ :— nR{n,A). 
Moreover (see |ENOO| Chapter II, Lemma 3.4]): 



Hm J„x — X for any x H, 

n — >oo 

Hm AnX — Ax for any x € D(A). 



(3.31) 



By Proposition II. 2 11 ASn{t)x = Sn{t)Ax for every n sufficiently large and for 
all X £ D{A). So, by Propositions 11.21] and and the closedness of A we 
can write 



AnW^it) = An Sn{t- T)1r{T)dW{T) 
JO 

^Jnf ASnit - T)^{T)dW{T) = J„ / Sn{t - T)A^{T)dW{T) 

Jo Uo 

Analogously, 

AnW^it) - J„ [ / Sit - T)A^iT)dWiT) 



By (13301) we have 
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NnAt) = \Jn / [Sn{t - t) - S{t - t)] A^f (T)dM^(T) |ff 



< 



[Snit - r) - S{t - T)]A^{T)dW{T)\H 



Since from assumptions A'l' S M^{Q,T; L2), then the term apearing above, 
[Sn{t — r) — S'(t — t)]A'1'{t) may be treated hke the difference defined by 

Hence, from p.30|) and (j3.28p . for the first term of the right hand side of 
([3:291) we have 

hm sup E{N^ ^{t)) ^ 0. 

For the second term of (|3.29p . that is 2(^)1 '^'^ can follow the same steps 
as above for proving (|3.28p . 

N^At)^\A„W''{t)~-AW'^{t)\H 

EE An [ S{t - T)<Ir{T)dW{T) -A[ S {t - t)>I' {T)dW (t) 
Jo JO 

= f [An- A]Sit - TpiT)dWiT) . 
JO H 

From assumptions, W,A^ G AP{0,T;L^). Because An,S{t), t > are 
bounded, then A„S{t- ■)^{-) £ 7V2(0,T;X§), too. 

Analogously, AS {t ~ {■) = Sit--)AW{-) e JV^iO,T; L°). 

Let us note that the set of all Hilbert-Schmidt operators acting from one 
separable Hilbert space into another one, equipped with the operator norm de- 
fined on page [20] is a separable Hilbert space. Particularly, sum of two Hilbert- 
Schmidt operators is a Hilbert-Schmidt operator, see e.g. |Ba81| . Therefore, 
we can deduce that the operator (A„ - A) S(t - ■)^{-) G A/'2(0, T; for 
any t G [0, T]. Hence, the term [An — A]S{t — t)^{t) may be treated like the 
difference defined by (|3.22|) . So, we obtain 



H 




^|K-A]5(t-rMr).g,|^ 



dr 



[[An- A]S{t-T)^{T)g,[^^dT\ <+cx). 



for any t G [0,T]. 

By the convergence p.3ip . for any fixed j, 

/ [[An-A]S{t-T)^{T)g,[]jdT^Q for 

JO 
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Summing up our considerations, we have 



Hm sup EiNl^{t))^0. 



te[o,T] 



So, we can deduce that 



hm sup E\AnW^{t) - AW'^it)\jf = 0, 



and then ([X^ holds. □ 

These considerations give rise to the following result. 

Theorem 3.16. Suppose that assumptions of Lemma \3.15\ hold. Then the 
equation (3. 1\) has a strong solution. Precisely, the convolution defined 
by 113.5]) is the strong solution to i3.1\) . 

Proof. Since Proposition l3.6l and Lemma [3.151 hold . we have to show only the 
condition (|3.2p . Let us note that by Proposition 13.71 the convolution W'^{t) 
has integrable trajectories. Because the closed unbounded linear operator A 
becomes bounded on see e.g. |We80) , we obtain that AW^l-) £ 

L^{[0, T]; H), P-a.s. Next, properties of convolution provide integrability of the 
function a(T — ■)AW'^ what finishes the proof. □ 

3.3 Fractional Volterra equations 

Assume, as previously, that H \s a. separable Hilbert space with a norm \ ■ \h 
and A\.s& closed linear operator with dense domain D{A) C H equipped with 
the graph norm | • |_d(a)- The purpose of this section is to study the existence 
of strong solutions for a class of stochastic Volterra equations of the form 



appropriate stochastic processes. There are several situations that can be 
modeled by stochastic Volterra equations, see e.g. [H096[ Section 3.4 ] and 
references therein. A similar equation was studied in [BT03| . too. Here we are 
interested in the study of strong solutions when equation (|3.32p is driven by 
a cylindrical Wiener process W . We give sufficient conditions for stochastic 
convolution to be a strong solution to ()3.32p . 

The equation ()3.32|) is a stochastic version of the deterministic Volterra 
equation 




where aa{t) : 



r{a) 



, a > 0, r{a) is the gamma function and W, ^Z' are 




(3.33) 



where / is an iJ-valued function. 
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In the case when aa{t) is a completely positive function, sufficient con- 
ditions for existence of strong solutions for p.32p may be obtained like in 
section 13. 2[ that is, using a method which involves the use of a resolvent 
family associated to the deterministic version of equation (j3.32p . 

However, there are two kinds of problems that arise when we study (j3.32p . 
On the one hand, the kernels are a-regular and ^-sectorial but not 



completely positive functions for a > 1, so e.g. the results in |KL07cj can 
not be used directly for a > 1. On the other hand, for a £ (0, 1), we have a 
singularity of the kernel int = 0. This fact strongly suggests the use of a-times 
resolvent families associated to equation (|3.33p . These new tools appeared in 
[BaOl] as well as their relationship with fractional derivatives. For convenience 
of the reader, we provide below the main results on a-times resolvent families 
to be used in this paper. 

Our second main ingredient to obtain strong solutions of p.32p relies on 
approximation of a-times resolvent families. This kind of result was very re- 
cently formulated by Li and Zheng ^LZ04) . It enables us to prove a key result 
on convergence of a-times resolvent families (see Theorem 13 . 201 below*) . Then 
we can follow the methods employed in |KL07c| to obtain existence of solu- 
tions - particularly strong - for the stochastic equation (|3.32p . 

3.3.1 Convergence of a-times resolvent families 

In this section we formulate the main deterministic results on convergence of 
resolvents. 

By Sa{t)., i > 0, we denote the family of a-times resolvent families corre- 
sponding to the Volterra equation p.33p . if it exists, and defined analogously 
like resolvent family, see Definition ll.il 

Definition 3.17. (see JBaOlf ) 

A family {Sa{t))t>o of bounded linear operators in a Banach space B is called 
a-times resolvent family for hS.SS]) if the following conditions are satis- 



1. Sa{t) is strongly continuous on R+ and Sa{0) — I; 

2. Sa{t) commutes with the operator A, that is, Sa{t){D{A)) C D{A) and 
ASait)x = Sait)Ax for all x e D{A) and t > 0; 

3. the following resolvent equation holds 



fied: 




(3.34) 



for all X £ D{A), t > 0. 



Necessary and sufficient conditions for existence of the a-times resolvent 
family have been studied in [Ba01| . Observe that the a-times resolvent fam- 
ily corresponds to a Co-semigroup in case a = 1 and a cosine family in case 
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a — 2. In consequence, when 1 < a < 2 such resolvent famihes interpolate 
Co-semigroups and cosine functions. In particular, for A = A, the integrodif- 
ferential equation corresponding to such resolvent family interpolates the heat 
equation and the wave equation, see |Fu89| or |SW89| . 

Definition 3.18. An a-times resolvent family {Sa(t))t>Q is called exponen- 
tially bounded if there are constants M > 1 and uj > such that 

\\Sa{t)\\ < Me'^\ t > 0. (3.35) 

If there is the a-times resolvent family {Sa{t))t>a for A and satisfying 
dSSni), we write A e C"(M,w). Also, set C"(w) 7= Um>iC"(M,w) and 
:= U^>oC"(t^). 

Remark It was proved by Bazhlekova [BaOli Theorem 2.6] that if A E C" for 
some a > 2, then A is bounded. 

The following subordination principle is very important in the theory of 
a-times resolvent families (see |Ba01[ Theorem 3.1]). 

Theorem 3.19. Let < a < /3 < 2, 7 = a/(3,uj > 0. If A e C^{lu) then 
A e C°'{uj^^'^) and the following representation holds 

/>oo 

Sa{t)x = / ipt,j{s)Sfj{s)xds, t > 0, (3.36) 
"'0 

where (pt^^{s) := t~^$^{st^^) and 'I'-^iz) is the Wright function defined as 

^>^(z) := y , ^, ^"""^ ^ -, < 7 < 1. (3.37) 



Remarks 1. We recall that the Laplace transform of the Wright function 
corresponds to E^{—z) where denotes the Mittag-Leffler function. In par- 
ticular, ^-yiz) is a probability density function. 

2. Also we recall from |Ba011 (2.9)] that the continuity in t > of the 
Mittag-LefHer function together with the asymptotic behavior of it, imply 
that for u! > there exists a constant C > such that 

L;„(wr) < Ce"'^°*, t>0, ae(0,2). (3.38) 

As we have already written, in this paper the results concerning conver- 
gence of a-times resolvent families in a Banach space B will play the key role. 
Using a very recent result due to Li and Zheng |LZ04| we are able to prove 
the following theorem. 
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Theorem 3.20. Let A be the generator of a Co-semigroup {T{t))t>o o- 
Banach space B such that 

||T(t)|| < Me'^*, t>0. (3.39) 

Then, for each < a < 1 we have A £ C°'{M,uj^/"'). Moreover, there exist 
hounded operators An and a -times resolvent families Sa,n{t) for An satisfying 
\\Sa.n{t)\\ < A/Ce(2")'^°*, for allt>0, ne N, and 

Sa,n{t)x — > Sa{t)x aS n +00 (3.40) 

for all X £ B, t > 0. Moreover, the convergence is uniform in t on every 
compact subset o/M_|_. 

Proof. Since A is the generator of a Cq semigroup satisfying (|3.39p . we have 
A G C^{uj). Hence, the first assertion follows directly from Theorem l3.191 that 
is, for each < a < 1 there is an a-times resolvent family {Sa{t))t>o for A 
given by 



Sa,it)x = / ipt,a{s)T{s)xds, t > 0. (3.41) 

Since A generates a Co-semigroup, the resolvent set p{A) of A contains 
the ray [w, oo) and 

\\R{X,A)''\\<j^-^^ ioYX>w, fceN. 

Define 

An ■■= nAR{n, A) = n^R{n, A) -nl, n> w, (3.42) 

the Yosida approximation of A. 
Then 

OO 2k-i-k 

I \e'^^ 1 1 = e-"* I |e"'«("'^)* 1 1 < e^"* ^ 1^^"' I 

fc=0 

Hence, for n> 2w we obtain 

||e^"*|| < Me^"'*. (3.43) 



Next, since each An is bounded, it follows also from Theorem 13.191 that for 
each < a < 1 there exists an a-times resolvent family {Sa,n{t))t>Q for An 
given as 



/•oo 

So..n{t) = / ¥'t,a(s)e"^"ds, t > 0. 
Jo 



(3.44) 
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By p.43p and remark 1, page [SOI it follows that 

/•oc /-oo 

Jo Jo 

= M ^a{T)e^'^'''''dT ^ MEa{2uje), t>0. 
Jo 

This together with remark 2, page [501 gives 

\\SaAt)\\ < AfCe(2")'^°*, t > 0. (3.45) 

Now, we recall the fact that i?(A,A„)x R{X,A)x as n ^ cx) for all A 
sufficiently large (see e.g. |Pa83[ Lemma 7.3]), so we can conclude from |LZ04[ 
Theorem 4.2] that 

Sa,nit)x ^ Sa{t)x SlS n ^ +00 (3.46) 

for all a; G B, uniformly for t on every compact subset of M_|_. □ 

An analogous result can be proved in the case when A is the generator of 
a strongly continuous cosine family. 

Theorem 3.21. Let A be the generator of a Co-cosine family {T{t))t>o in 
a Banach space B. Then, for each < a < 2 we have A g C"(M, w^/"). 
Moreover, there exist bounded operators An and a-times resolvent families 
Sa,n{t) for An Satisfying \ \So,,n{t)\\ < MCe^^'^^''"* , for allt>0, n e N, and 

Sa.7i{t)x Sait)x as n — > +00 

for all X € B, t > 0. Moreover, the convergence is uniform in t on every 
compact subset 0/IR.+ . 

Let us note that formulae (|3.41|) and l|3.44p still hold when A is the Cq- 
cosine family and < a < 2. 

In the following, we denote by S0{uj) the open sector with vertex w e M 
and opening angle 29 in the complex plane which is symmetric with respect 
to the real positive axis, i.e. 

Seiuj) := {A e C : \arg{\ - lu)\ < 9}. 

We recall from [BaOL, Definition 2.13] that an a-times resolvent family 
Sa(t) is called analytic if Sa{t) admits an analytic extension to a sector 
Sdg for some Oq G (0, 7r/2]. An a-times analytic resolvent family is said to 
be of analyticity type {9q,luq) if for each 9 < 9o and to > ujq there is 
M = M (9, oj) such that 
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The set of all operators A ^ C" generating a-times analytic resolvent fam- 
ilies Sa{t) of type (6'o,wo) is denoted by A°'{9o,uJa)- In addition, denote 
^"(^o) := UM"(^o,t^o);^o G ^" := U{^"(^o);eo e (0,^/2]}. For 

a = 1 we obtain the set of all generators of analytic semigroups. 

Remark We note that the spatial regularity condition TZ{Sa{t)) C D{A) for 
alH > is satisfied by a-times resolvent families whose generator A belongs 
to the set ^"(^Oj'^o) where < a < 2 (see jBaOll Proposition 2.15]). In 
particular, setting = we have that A S A"{9o,0) if and only if —A is 
a positive operator with spectral angle less or equal to tt — a{TT/2 + 6). Note 
that such condition is also equivalent to the following 

^a(V2+e) and ||A(A/-A)-i|| < Af, A e r„(,/2+e). (3.47) 

The above considerations give us the following remarkable corollary. 

Corollary 3.22. Suppose A generates an analytic semigroup of angle 7r/2 and 
a G (0, 1). Then A generates an a-times analytic resolvent family. 

Proof. Since A generates an analytic semigroup of angle 7r/2 we have 

\\X{XI~A)-^\\<M, XeS^-,. 

Then the condition (|3.47p (see also |Ba01[ Corollary 2.16]) implies A e 
^"(min{^^7r, i7r}, 0), a £ (0,2), that is A generates an a-times analytic 
resolvent family. □ 

In the sequel we will use the following assumptions concerning Volterra 
equations: 

(Al) A is the generator of Co-semigroup in H and a € (0, 1); or 
(A2) A is the generator of a strongly continuous cosine family in H and 
a e (0,2). 

Observe that (A2) implies (Al) but not vice versa. 
3.3.2 Strong solution 

As previously H and U are separable Hilbert spaces and W is a cylindri- 
cal Wiener process defined on a stochastic basis {f^,T, {J-')t>o, P), with the 
positive symmetric covariance operator Q e L{U), TrQ = +oo. The spaces 
Uq, = L2{Uq, H) and A/'(0, T; L^) are the same like in previous sections. 

For the reader's convenience we formulate definitions of solutions to the 
equation (|3.32p . We define solutions to the equation (|3.32p analogously like in 
section 13.11 
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Definition 3.23. Assume that (PA) hold. An H -valued predictable process 
X{t), t £ [0,T], is said to be a strong solution to \3.S2]) . if X has a 
version such that P{X{t) G D{A)) = 1, for almost all t G [0,T]; for any 
tG[0,T] 

t 

\aa{t- t) AX {t)\h dr < P-a.s., a > 0, (3.48) 
and for any t G [0, T] the equation US. 32\) holds P — a.s. 

Definition 3.24. Let (PA) hold. An H -valued predictable process X{t), t G 
[0,r], is said to be a weak solution to \3.32]) . if P{Jq \aait - T)X{T)\HdT < 
-\-oo) = 1, a > 0, and if for all ^ G D{A*) and all t G [0,r] the following 
equation holds 

{X{t), Oh = (^(0), OH + {f a^{t- t)X{t) dr, A*Oh 

Jq 

+ ( / <f(r)dVK(r),e)H, P-a.s. 
Jo 

Definition 3.25. Assume that X{0) is J-q -measurable random variable. An 
H-valued predictable process X(t), t G [0,r], is said to be a mild solution 
to the stochastic Volterra equation h3.32]) . if E(/g ||S'a(t - r)if'(T)||^„ dr) < 
+CXD, a > 0, for t <T and, for arbitrary t G [0, T], 

X{t) = So,{t)X{0) + [ Sait-T)W{T)dW{T), P-a.s. (3.49) 

where Sa{t) is the a-times resolvent family. 
We define the stochastic convolution 

W^it):^ f 5a(t-T)tf'(r)rfW^(r), (3.50) 
Jo 

where G J\f^{0,T; L2). Because a-times resolvent families Sa{t), t>0, are 
bounded, then S^it - G A/'2(0, T; L^), too. 

Analogously like in section [^TTl we can formulate the following results. 

Proposition 3.26. Assume that Sa{t),t > 0, are the resolvent operators to 
iTOi)) . Then, for any process if G 7V2(0,T;i§), the convolution W^{t), t > 
0, a > 0, given by 113. 5\) has a predictable version. Additionally, the process 
W^{t), t > 0, a > 0, has square integrable trajectories. 

Under some conditions every mild solution to (|3.32p is a weak solution to 
(13321) . 
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Proposition 3.27. If ^ e Af^{0,T; L9,), then the stochastic convolution 
fulfills the equation 

«W,Oh= f\a^{t-T)W^{T),A*OH+ f\^Mr)dW{T))H, (3.51) 
Jo Jo 

a>0 for any t G [0,T] and ^ € D{A*). 

Immediately from the equation (j3.5ip we deduce the following result. 

Corollary 3.28. If A is a bounded operator and <F e 7V^(0,T;i^), then the 
following equality holds 

W^it)^ f a^{t~T)AW^{T)dT+ [ ^{T)dW{T), (3.52) 
Jo Jo 

forte [0,T], a>0. 

Remark The formula (|3.52p says that the convolution W^{t), t > 0, a > 0, 
is a strong solution to (|3.32p if the operator A is bounded. 

We can formulate following result which plays a key role in this subsection. 

Lemma 3.29. Let assumptions (VA) he satisfied. Suppose that (Al) or (A2) 
holds. If^ and A^ belong to Af^{0,T; L^) and in addition ^{t){Ua) C D{A), 
P-a.s., then the following equality holds 

Wl{t)^ [ a„(t-r)AW^f (T)dT+ / ^{T)dW{T), P-a.s. (3.53) 
Jo Jo 

Remark Let us emphasize that in (Al), a e (0, 1) and in (A2), a £ (0, 2). 

Although the proof is analogous to that given in section [231 we formulate 
it for the reader's convenience. 

Proof. Because formula (|3.52p holds for any bounded operator, then it holds 
for the Yosida approximation An of the operator A, too, that is, 

f aa{t~T)A^WljT)dT+ f^{T)dW{T), 

Jo Jo 

where ^ 

WZnit) ■■= I S^At-rmr)dW{r) 
Jo 

and ^ 

AnW:Zjt)=An f S^,n{t-Tp{T)dWiT). 

Jo 

By assumption •f' e M^{0,T; L2). Because the operators Sa,n{t) are deter- 
ministic and bounded for any t G [0, T], a > 0, n £ N, then the operators 
Sa,n{t — ■)'!'{■) belong to A/'^(0, T; L2), too. In consequence, the difference 
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^cAt - •) - S^At - -Pi-) - S^it - ■)*■(•) (3-54) 
belongs to 7V^(0, T; ij) for any t e [0, T], a > and n e N. This means that 



E 



(3.55) 



for any t G [0,T]. 

The cylindrical Wiener process W{t), t>0, can be expanded in the series 



+ 00 



(3.56) 



where {gj} is an orthonormal basis of Uq and f3j{t) are independent real 
Wiener processes. From p.56p we have 



/ ^c.,n{t-r)dW{T)^J2 '^aAt-T)gjdPj{T). (3.57) 
Jo Jo 



In consequence, from p.55p 
„t / +00 



E 



< +00 



(3.58) 



for any t £ [0,T]. Next, from (|3.57p . properties of stochastic integral and 
(|338ll we obtain for any t £ [0, T], 



E 



<?a,„(i-T) dW{T) 



= E 



+00 „t 

^ / <Z'a.„(i-T)g,d/3,(T) 
,=1 -^0 

+ 00 „T 

\<^aAT-r)g,\Ur 

.7=1 



E 



5Z / |^>a,„(i-T)5,|^dT 



< E 



< 



< +00. 



By Theorem 13.201 or Theorem I3.21[ the convergence of a-times resolvent 
families is uniform in t on the interval [0, T]. So, for any fixed a and j, 

\[SaAT-T)~ Sa{T-T)]^{T)gj\'jjdT — >0 for n->+oo. (3.59) 

Then, using ()3.58|) and (|3.59|) we can write 

2 



sup 

te[o,T] 



$c.At-^)dW{T) 



H 



= sup E 

te[o,T] 







3.3 Fractional Volterra equations 

t 2 
[Sa.nit -T)-Sa{t- T)]<F{T)dW{T) 

H 
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< E 



TOO „T 

E / \[S^,n{T-r)-SUT~rmr)9j\ldr 
,=1 -^0 



as n — > +00 for any fixed a > 0. 

Hence, by the Lebesgue dominated convergence theorem we obtained 



liin sup E\wl^{t)-W:^{t)\^=0. 



(3.60) 



By Proposition EH (i) G D{A)) = 1. 

For any n G N, a > 0, t > 0, we have 



where 



Then 



(3.61) 



Let us study the term Nn.i{t)- Note that, either in cases (Al) or (Al) the 
unbounded operator A generates a semigroup. Then we have from the Yosida 
approximation the following properties: 

AnX = JnAx for any x e D{A), sup 1 1 J„| | < oo, (3.62) 

n 

where AnX = nAR(n,A)x — AJnX for any x Cz H with J„ := nR{n,A). 
Moreover (see [ENOOl Chapter II, Lemma 3.4]): 

lim n R{n, A)x — x for any x G iJ, (3.63) 

n — >oo 

lim AjiX = Ax for any x G D{A). (3.64) 

n — *oo 

Note that ASa,n{t)x = Sa,n{t)Ax for all x G D{A), since e*"^" commutes 
with A and A is closed (see p.44p ). So, by Proposition 13.141 and again the 
closedness of A we can write 

AnWo,,„{t) = A„ I 5a,„(t-T)lf(T)dW^(T) 



Sc.,r^{t~T)A^{T)dW{T) 



nR{n, A) 
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Analogously, 

AnWait) = nR{n, A) f Sc,{t - T)A^{T)dW{T) 
Uo 

By (|3.63p we have 

iV„.l(i) = I J„ / [So. At -T)~Sa.{t~ T)]A^{T)dW{T)\H 

Jo 

< I / [Sa.nit - r) - So.it - T)]A>F{T)dW{T)\H . 

Jo 

From assumptions, A^ e Af'^{0,T; L^). Then the difference [6*0, „(t ~ t) ~ 
Sa{t — t)]A'I'{t) may be estimated exactly like the difference 'Pa,n defined by 
([334)l . 

Hence, from p.63p and ()3.60p for the first term of the right hand side of 
(|3.6ip we obtain 

lim sup E{N^ -> 0. 

For the second and third terms of (j3.61|) we can follow the same steps as above 
for proving ()3.60|) . We have to use the properties of Yosida approximation, 
particularly the convergence (|3.64p . So, we can deduce that 

hrn sup m^W'i^it) - AW^it)\j, = 0, 

what gives p.Sp . □ 

Now, we are able to formulate the main result of this section. 

Theorem 3.30. Suppose that assumptions of Lemma \3. 29\ hold. Then the 
equation IjS. 32]) has a strong solution. Precisely, the convolution defined 
by i3.50\) is the strong solution to h3.3^) . 

Proof. We have to show only the condition (|3.48p . The convolution W^^t) has 
integrable trajectories (see section [3?T|) . that is, W^(-) e L^{[Q,T]]H), P-a.s. 
The closed linear unbounded operator A becomes bounded on {D{A), \ ■ \d{A))i 
see ^We80. Chapter 5]. So, we obtain AW^[-) € L^[[Q,T]-H), P-a.s. Hence, 
the function aa{T — t)AW'^{t) is integrable with respect to r, what finishes 
the proof. □ 

The following result is an immediate consequence of Corollary 13.221 and 
Theorem ESOl 

Corollary 3.31. Assume that (VA) hold, A generates an analytic semigroup 
of angle 7r/2 and a S (0,1). If ^ and A^ belong to 7V^(0, T;L^) and in 
addition •f'(i)([/o) C D[A), P — a.s., then the equation i3.38\) has a strong 
solution. 
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3.4 Examples 

In this short section we give several examples fulfilling conditions of theorems 
providing existence of strong solutions. The class of such equations depends 
on where the operator A is defined, in particular, the domain of A depends 
on each considered problem, and also depends on the properties of the kernel 
function a{t), t >0. 

Let G be a bounded domain in M" with smooth boundary dG. Consider 
the differential operator of order 2m: 

A{x,D)= ba{x)D" (3.65) 

\a\<2m 

where the coefiicients ha{x) are sufficiently smooth complex- valued functions 
of X in G. The operator A{x^D) is called strongly elliptic if there exists 
a constant c > such that 

i?e(-l)™ > cl^P" 

\a\=2m 

for all a; e G and ^ e M". 

Let A(x^ D) be a given strongly elliptic operator on a bounded domain 
G C M" and set D{A) = H^'^iG) n H^{G). For every u e D{A) define 

Au = A{x,D)u. 

Then the operator —A is the infinitesimal generator of an analytic semigroup 
of operators on fi" = L'^{G) (cf. |Pa83[ Theorem 7.2.7]). We note that if the 
operator A has constant coefficients, the result remains true for the domain 

G = M". 

The next example is the Laplacian 

2—1 * 

since —A is clearly strongly elliptic. It follows that Au on D{A) — H'^{G) n 
Hq{G) is the infinitesimal generator of an analytic semigroup on L^{G). 

In particular, by [Pr93, Corollary 2.4] the operator A given by p.65p gen- 
erates an analytic resolvent S{t) whenever a G G(0, oo)nL^(0, 1) is completely 
monotonic. 

This example fits in our results iia{t) is also completely positive. For exam- 
ple: a{t) ~ t"^^ / r{a) is both, completely positive and completely monotonic 
for < a < 1 (but not for a > 1). 

Another class of examples is provided by the following: suppose a £ 
L[q^(R+) is of subexponential growth and 7r/2-sectorial, and let A gener- 
ate a bounded analytic Gp-semigroup in a complex Hilbert space H. Then 
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it follows from [Pr93[ Corollary 3.1] that the Volterra equation of scalar 
type u = a * Au + / is parabolic. If, in addition, a{t) is fc-regular for all 
fc > 1 we obtain from |Pr93i Theorem 3.1 ] the existence of a resolvent 
S e C'="i((0,oo),S(i?)) such that 11(3(1)) c D{A) for i > (see [P793l 
p.82 (f)]). 
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Stochastic Volterra equations in spaces of 
distributions 



In this chapter we study two classes of hnear Volterra equations driven by 
spatially homogeneous Wiener process. We consider existence of solutions to 
these equations in the space of tempered distributions and then derive condi- 
tions under which the solutions are function-valued or even continuous. The 
conditions obtained are expressed in terms of spectral measure and the space 
correlation of the noise process, as well. Moreover, we give description of 
asymptotic properties of solutions. 

The chapter is organized as follows. In section H?T] we introduce generalized 
and classical homogeneous Gaussian random fields basing on |GV64| . [AdSlj 
and [PZ97| . We recall some facts which connect the generalized random fields 
with their space corralations and spectral measures. Moreover, we recall some 
results used in the proofs of the main theorems. Section [42l originates from 
|KZ00aj . Here we study regularity of solutions to the equation (|4.ip and pro- 
vide some applications of these results. Section |473] is a natural continuation of 
section [4.21 In this section we give necessary and sufficient conditions for the 
existence of a limit measure to the stochastic equation under consideration. 
Results of section 14.31 come from [Ka03 . In section 14.41 we study an integro- 
differential stochastic equation with infinite delay. We provide necessary and 
sufficient conditions under which weak solution to that equation takes values 
in a Sobolev space. Section l44l originates from |KL07aj . 

4.1 Generalized and classical homogeneous 
Gaussian random fields 

We start from recalling several concepts needed in this chapter. 

Let SiR'^) , SciR'') , denote respectively the spaces of all infinitely differ- 
entiable rapidly decreasing real and complex functions on Mf^ and S"(R'*), 
iS'^(R'^) denote the spaces of real and complex, tempered distributions. The 
value of a distribution ^ e S'^{M.'^) on a test function ip will be written as 
(^, '(/')■ For ?A G SiW^) we set ^(s){0) = ipi-O), 9 S W^. Denote by 5'(,)(R'*) 
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the space of all ^ G S{R'^) such that i/j = V'(s), and by 5*1^-1 (M'') the space of 
aU C e S'{R'^) such that {C,Tp) = (t^is)) for every 4> e S{R'^). 

We define the derivative ^ of the distribution ^ e 5'(M'') by the formula 
(l tp) = -(^, for if € S'(M''), see [GS64] . 

In the chapter we denote by ^ the Fourier transform both on S'c(K'*), and 
on S'^iW^). In particular, 

and for the inverse Fourier transform 

Moreover, if ^ G S'^(R'^), 

for all V' G 5c (M''). Let us note that JF transforms the space of tempered 
distributions S'{R'^) into S[^^(R'^). 

For any h gR'^, ip e S{R'^), C e S"(R'^), the translations r,i?/>, r;;^ are 
defined by the formulas 

By ;B(5'(R'^)) and S(5^(R'^)) we denote the smallest a-algebras of subsets 
of 5"(R'') and ^^(R''), respectively, such that for any test function ip the 
mapping — > f) is measurable. 

The below notions of generalized random fields, their space correlations 
and spectral measures are recalled directly from |GV64| . 

Let J^, P) be a complete probability space. Any measurable mapping 
Y : fi ^ S'iW^) is called a generalized random field. A generalized 
random field Y is called Gaussian if {Y,ip) is a Gaussian random variable 
for any (p € 5'(R''). The definition implies that for any functions ipi, . . . , (pn G 
S'(IR.'^) the random vector {{Y, ipi), . . . ,{Y, (/?„)) is also Gaussian. One says that 
a generalized random field Y is homogeneous or stationary if for all h G R'^, 
the translation rj^iY) of Y has the same probability law as Y. 

A distribution F on the space 5(R'') is called positive-definite if 
{r,(p -k ip(^g-j) > for every ip G S{R'^), where ip -k ip^g^ denotes the convolu- 
tion of the functions ip and (p(s)- 

If F is a homogeneous, Gaussian random field then for each i/j G S{R'^), 
(Y, t/j) is a Gaussian random variable and the bilinear functional q : S(R'^) x 
S{R'^) R defined by the formula. 



=IE((r,^)(r,^)) for ^e5(M'^), 
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is continuous and positive definite. Since q{Lp,^) = q{Tii^,Tfiip) for all ip, 
ip € S'(M'^), h G W^, there exists, see e.g. (GV641 Chapter II], a unique positive- 
definite distribution F € S' {W^) such that for all (p, il> € S{M.'^) , one has 

The distribution F is called the space correlation of the field Y. By 
Bochner-Schwartz theorem the positive-definite distribution F is the inverse 
Fourier transform of a unique positive, symmetric, tempered measure fi on 
R'' : 7^ = JF^^(/i). The measure is called the spectral measure of F and 
of the field Y. 

Summing up a generalized homogeneous Gaussian random field Y is char- 
acterized by the following properties: 

1. for any ip € S'(R'^), {Y,^fj) is a real-valued Gaussian random variable, 

2. there exists a positve-definite distribution 7^ e S"(R'') such that for all 

E((y,^)(y,7^)) = (r,¥.*^(,)), 

3. the distribution F is the inverse Fourier transform of a positive and sym- 
metric tempered measure fj, on Mf^, that is, such that 

(1 + |A|)'' fi{dX) < +00, for some r < . 

Let Y : Q ^ S'{W^) be a generalized random field. When the values of Y 
are functions, with probability 1, then Y is called a classical reindom field 
or shortly random field. In this case, by Fubini's theorem, for any 9 € 
the function Y(9) is well-defined and 

(r,<^)= / Y{e)^{e)de, 

for any </? G S'(K''). Thus any random field Y may be identified with the 
family of random variables {Yg^g^^d parametrized by G . In particular 
a homogeneous (stationary), Gaussian random field is a family of Gaussian 
random variables Y{9), 9 e M"*, with Gaussian laws invariant with respect 
to all translations. That is, for any 9i,. . . ,9n G R'' and h G Mf^, the law of 
{Y{9 -f /i), . . . , Y{9„ + h)) does not depend on h e R^. 

For the sake of completness we sketch now the proof of the following result, 
(see also [PZ97] ). 

Proposition 4.1. A generalized, homogeneous, Gaussian random field Y is 
classical if and only if the space correlation FofYisa bounded function and 
if and only if the spectral measure ^ of Y is finite. 

Proof. First, let us prove that if a positive definite distribution r' is a bounded 
function then it is continuous and its spectral measure /z is finite. It is enough 
to show only that the spectral measure /i of the distribution F is finite. 



Jm 
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Let pt(-) denote the normal density with the Fourier transform e , 
t > 0, A e R''. Define measures /it, t > 0, by the formula 

J B 

Since the measure /i is tempered, the measures /it are finite. 

The Fourier transform T{^t) of /if, for any t > 0, is a continuous function and 



But 



and 



^t{e) = T-\^Hm = {^*ptm, i 
rt{0)= [ r{0~T^)pt{r,)dTi 



e 



sup |r(C)| 



sup inoi 



In particular 

Letting 1 1 one obtains that 



n{dX) < 



sup \ric)\ 



sup |r(c)| 



< +00 , 



so the measure fi is finite as required. 

Let now F be a classical, homogeneous, Gaussian random field. It means 
that the field Y is function-valued. Moreover E{Y{ei)Y{92)) = r{0i - 62), 
for 9i, 02 £ M'', the space correlation F is positive-definite and \r{9i — 62)1 < 

r(o) < +00. 

Let now fi be the finite spectral measure of a homogeneous Gaussian 
random field Y. Then F is a positive definite continuous function. By Kol- 
mogorov's existence theorem, there exists a family Y{d), 9 G M'^, such that: 



E(y(0i)y(02)) = m 



h,92 e 



From the continuity of F, it follows that the family Y, is stochastically con- 
tinuous and therefore has a measurable version. Since the laws of the random 
fields Y, Y coincide, the result follows. □ 

We finish the section recalling a continuity criterium which will be used in 
the proof of the continuity results, see f [Ad81| . Th. 3.4.3). 

Proposition 4.2. Let Y{9), 9 e R'', be a homogeneous, Gaussian random 
field with the spectral measure fj,. If, for some e > 0, 



(ln(l + |A|))i+^/i(dA) <+oo. 



then Y has a version with almost surely continuous sample functions. 
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4.2 Regularity of solutions to stochastic Volterra 
equations 

This section is concerned with the following stochastic Volterra equation 

x{t,e)^Xo{0)+ [ b{t~T)Ax{T,e)dT + w{t,e), (4.1) 

Jo 

where t eR+, 9 £ R'^, Xq e 5'(M''), b e Ll^{M.+) and is a spatially ho- 
mogeneous Wiener process which takes values in the space of real, tempered 
distributions S"(M''). The class of operators A covered in the present chap- 
ter contains in particular the Laplace operator A and its fractional powers 
-{-A)^/^ p e {0,2]. ^ 

The equation l|4.ip is a generalization of stochastic heat and wave equations 
studied by man y au thors, s ee e.g. |DF98| . [KZ01| . [KZOOj . |MMOO| . |MS99| . 
|Mu97| . |PZOO] and [Wa86| and references therein. In the context of infinite 
particle systems stochastic heat equation of a similar type has been investi- 
gated by Bojdecki with Jakubowski [B J90| . [BJ97j and |BJ99j and by Dawson 
with Gorostiza in [DG90| . 

As we have already said, our aim is to obtain conditions under which 
solutions to the stochastic Volterra equation (I4.1|l are function-valued and 
even continuous with respect to the space variable. In the chapter we treat 
the case of general dimension and the correlated, spatially homogeneous noise 
Wr of the general form. 

4.2.1 Stochastic integration 

In this chapter we will integrate operator-valued functions Tl{t), t > 0, with 
respect to a Wiener process W. The operators Tl{t), t > 0, will be non-random 
and will act from some linear subspaces of S"(R'^) into S"(R''). We shall assume 
that W{t), t > 0, is a continuous process with independent increments taking 
values in S"(R'^). The process W is space homogeneous in the sense that, 
for each t > 0, random variables W{t) are stationary, Gaussian, generalized 
random fields. We denote by F the covariance of W{1) and the associated 
spectral measure by ^. To underline the fact that the probability law of W 
is determined by F we will write Wr- From now on we denote by g a scalar 
product on S'(R'^) given by the formula 

(j(^,^) = (r,^*^(,)), ^^^eSiR"). 

Let us present three examples of spatially homogeneous Wiener processes. 
Examples 1. Important examples of random fields are provided by symmetric 
a-stable distributions F{x) = e"'^' , where a s]0, 2]. For a = 1 and a = 2 the 

_ d + l 

densities of the spectral measures are given by the formulas ci(l-f-|a;p) ^ 
and C2e~l^l , where ci and C2 are appropriate constants. 
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2. Let q{'4>, ip) — (^{^A + to^) ^ ip,ip^, where A is the Laplace operator on 
and TO is a strictly positive constant. Then r' is a continuous function on 

E'*\{0} and ^{x) = {2Tr)-i {\x\^ + m^y\ The lawof V7(l) is the so-called 
Euclidean free field. 

3. Let q{ip,ip) = {'>p,ip). Then F is equal to the Dirac Nonfunction, its 



spectral density is the constant function (27r) 2 and is a white noise 
on L^([0,oo[xR''j. If B{t,x), t > and x e R"*, is a Brownian sheet on 



dx 

[0, oo[xR'^, then W can be defined by the formula 



, d'^B(t,x) 

The crucial role in the theory of stochastic integration with respect to 
Wr is played by the Hilbert space S'^ C S"(]R'') called the kernel or the 
reproducing kernel of Wp- Namely the space S'^ consists of all distributions 
^ e S'iW'-) for which there exists a constant C such that 



The norm in 5'^ is given by the formula 

Let us assume that we require that the stochastic integral should take values 
in a Hilbert space H continuously imbedded into S'{W^) . Let LnsiSg, H) be 
the space of Hilbert-Schmidt operators from S'^ into H . Assume that 7?,(t), 
t > is measurable Lhs{S'^, i?)-valued function such that 

/ l|7^('^)lli„s(S' i?)^^ < +°°' foraUt>0. 
Then the stochastic integral 

[ n{a)dWr i<7), t>0 
Jo 

can be defined in a standard way, see [It84| , |DZ92] or [PZ97| . The stochastic 
integral is an if-valued martingale for which 



E n{a)dWr{(Tij =0, t>0 



and 

nia)dWr{a) - E ( / ||7^(a)||i^^(5, ) , t > 0. 



E 



'0 



H 
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We will need a characterization of the space S'^. In the proposition below, 
L^j,j(M'', /i) denotes the subspace of L^(M'', /i; C) consisting of all functions u 
such that M(s)(6') u{-9) for 6 e R'^. 

Proposition 4.3. \PZ9'T\ Proposition 1.2] A distribution ^ belongs to S*^ if 
and only if ^ = ufl for some u <ELf^^{R'^,^i). Moreover, if ^ = ufl and rj — vjl, 
then 

The operators TZ{t), t > 0, of convolution type are of special interest 

TZ{t)£_ = r{t) * £_ , t>0, ee5'(M''), 

with r{t) e S'{M.'^) . The convolution operator is not, in general, defined for 
all ^ £ S'{M.'^) and for the stochastic integration it is important to know 
under what conditions on r(-) and ^ the convolution is well-defined. For many 
important cases the Fourier transform JFr(t)(A) , i > , A € R'', is continuous 
in both variables and, for any T > , 

sup sup \Tr{t){X)\ = Mt < +oo. (4.2) 

te[o,T] AeR'' 

If this is the case then the operators Tl{t) can be defined using Fourier trans- 
forms 

nit)^^T-'{Tv{t)H), 

for all ^ such that JF^ has a representation as a function. 

Now, we can characterize the stochastic convolution as follows. 

Theorem 4.4. Assume that the function J-r is continuous in both variables 
and satisfies condition \4--2^ . Then the stochastic convolution 

n*Wr{t)^ I n{t~ a)dWri(T) , t>0, 
Jo 

is a well-defined S'lW^) -valued stochastic process. For each t>0, TZ*Wr{t) 
is a Gaussian, stationary, generalized random field with the spectral measure 

l^tidX) = (^^ \Tr{a){X)\''da^ ^l{dX) , (4.3) 

and with the covariance 

Ft ^ f Y{a)*F*Y(s){,(T)dG . (4.4) 
Jo 
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Proof. Let p be an arbitrary continuous scalar product on S'(IR.'') , such that 
the embedding S'^iR'^) C ^^(R'') is Hilbert-Schmidt; here ^^(R'^) denotes 
a space of distributions on S{R'^) endowed with p. For more information on 
a family of Hilbert spaces of distributions we refer to [It83j . Note that for 

^^T'^iufi) with w G L^,)(M'^,^) . 
By (1121), T{r{t))T{C) is a measure 

^(r(t))(A)u(A)M(dA), 

belonging again to S'giR'^) . Moreover, for t G [0,T] , 

mmLiS'S') < sup sup \Tr{t){X)\ = Mt < +(X3. 

te[o,T] xeR'' 

Since the embedding S'^ C S'p is Hilbert-Schmidt, the stochastic integral, by 
the very definition, is an S'p -valued random variable. Denote 

Zt = n*Wr{t). 

Then we may write 



E((Zt,(^)(Zt,V^))=E(^(^y n{t-a)dWria),^j {^J n{t - u)dWriu),^ 
= [ {r{t-a)*(p,dWr{<7)) [ {r{t ~ u) * 'iP,dWr{u)) 



/ {r,{r{a) * ip) * {r{a) *^p)(s))da- 
Jo 

where (f, G S'(R'^) . This implies the formula (14. 4p of the theorem, from 
which (14. 3p easily folows. □ 



As an application define 



z;(A) = i(QA,A>- / {e^^^'y'' - l^dy) (4.5) 

In the formula (|4.5p . Q is a symmetric, non-negative definite matrix and v is 
a symmetric measure concentrated on \ {0} such that 

\y\'^iy{dy) < +oo, / liy{dy) < +00 . (4.6) 
y\<i J\y\>i 

This is the Levy-Khinchin exponent of an infinitely divisible symmetric law. 
From Theorem 14.41 we have the following proposition. 
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Proposition 4.5. Assume that 

J^r(t)(A) = e~*^(^\ t>0 

where v is the Levy-Khinchin exponent given by |^.5p and Then the 

conditions of Theorem \4-4\ o'^s satisfied. 

For more information on stochastic integral with values in the Schwartz 
space of tempered distributions S"(M'') we refere to Ito ( |It83j . [It84j ). Bo- 
jdecki with Jakubowski ([BJ89^-^BJ99^), Bojdecki with Gorostiza ( |BG86] ) 
and Peszat with Zabczyk (£Z97j-£Z00j|). 

4.2.2 Stochastic Volterra equation 

We finally pass to the linear, stochastic, Volterra equation in 5"(R'^) 

X{t)=Xo+ [ b{t-T)AX{T)dT + Wr{t), (4.7) 
Jo 

where Xq £ S"(R''), A is an operator given in the Fourier transform form 

HAOW^-vWmW, (4.8) 

V is a locally integrablc function and Wr is an 5'(]R'^)-valucd space homoge- 
neous Wiener process. 

Note that if v{X) = jAp, then A = A and if v{X) = |A|", a G (0,2), then 
A = — (— Z\)t is the fractional Laplacian. 

We shah assume the following HYPOTHESIS (H): 

For any 7 > 0, the unique solution s{-,"f) to the equation 

s{t)+'j [ b{t~T)s{T)dT ^1, t>0 (4.9) 

fullfils the following condition: for any T > , 

sup sup |s(f, 7)1 < +00. 
te[o,T] 7>0 

Comment Let us note that under assumption, that the function 5 is a locally 
integrable function, the solution s(-,7) of the equation (|4.9p is locally inte- 
grablc function and measurable with respect to both variables 7 > and t > 0. 

For some special cases the function s{t\ 7) may be found explicitely. 
Namely, we have (see e.g. |Pr93] ): 
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for b{t) = 1, s{t; 7) = e'^\ t>0, 7 > 0; 



(4.10) 



for b{t) = t, s{t; 7) = cos(/7i), i > 0, 7 > 0; 



(4.11) 



for b{t)=e-', s(t;7) = (l+7)-i 
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-(1+7)* 



t>0, 7 > 0. (4.12) 



We introduce now the so called resolvent family TZ{-) determined by the 
operator A and the function v . Namely, 

where, 

r{t) = T-'^s{Uv{-)), t>0. 

As in the deterministic case the solution to the stochastic Volterra equation 
((47)1 is of the form 

x{t) ^n{t)XQ + [ n{t^T)dWr{T), t>o. (4.13) 

Jo 

By Theorem 14. 4[ we can formulate the following result. 

Theorem 4.6. (\KZ00a[ Theorem 1]) Let Wp be a spatially homogeneous 
Wiener process and TZ(t), t > 0, the resolvent for the equation |^.7p . If Hy- 
pothesis (H) holds then the stochastic convolution 

n*Wr{t)^ [ n{t-a)dWr{a), t>0, 
Jo 

is a well-defined S^M."^) -valued process. For each t > the random variable 
7Z * Wr{t) is generalized, stationary random field on M.'^ with the spectral 
measure 



Ht{dX) 



t 

(s(cr, w(A)))^ da 



fi(dX). (4.14) 



By Propositions 14.11 and 14.21 we can conclude the below result. 

Theorem 4.7. ( \KZ00a\ Theorem 2]) A ssume that the Hypothesis (H) 
holds. Then the process TZ * Wr(t) is function-valued for all t > if and 
only if 

{s{a,v{X))f da] ^{dX) < +00 , t>0. 



\Jo 

If for some e > and all t > 



(ln(l + |A|))'+'' (s(cr, v{\))ydan{dX) < +oo , 
then, for each t > , TZ* Wr{t) is a sample continuous random field. 
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4.2.3 Continuity in terms of F 

In this subsection we provide sufficient conditions for continuity of the solu- 
tions in terms of the covariance kernel F of the Wiener process Wr rather 
than in terms of the spectral measure as we have done up to now. Analogi- 
cal conditions for existence of function-valued solutions can be derived in a 
similar way. 

Theorem 4.8. Assume that d > 2 , that F is a non-negative measure and 
A = -(-Z\)"/^ a e]0, 2] . If for some S>0, 



then for the cases _?_?[ ) and the solution to the stochastic 

Volterra equation |^.7| j has continuous version. 

The proof will be based on several lemmas. For any 7 G]0, 2] denote 
by pI the density of the 7-stable, rotationally invariant, density on the d- 
dimensional space. Thus, 

e-*l^r =.Fp7(A) . (4.15) 
Lemma 4.9. For arbitrary t > Q and arbitrary x G , 

p-l{x) ^t^'^/'^plixt-^'^). 
Proof. From (I4.15P we have 

Introducing a new variable y — t-i x, one has 

/: = / e'^^-y^pl{yr^)dy 
and the result follows. □ 

Lemma 4.10. There exists a constant c > such that for all "f < 2 , 



/ e-'pUx)dt<j-^, xe^" 



Proof. It is well-known, see e.g. Gorostiza and Wakolbinger |GW91l (p. 286)], 
that for some constant Ci > 0: 



Pji-)<Y^^, -GK"^- (4.16) 
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From Lemma 14.91 and the estimate (j4.16p we obtain: 



G2^) = 



+ CXD 



dt < 



+ 00 



+ —2. 



Cl 



< 



< 



+ 00 



Cl 



+ 00 



Cit~ 



t~ + \x\<^+'< 



dt < 



+ 00 



xt 

Cl 



<i+7 



dt 



i~ + |.-r|''+T 



dt 



i+7 



□ 



Lemma 4.11. 7/7 < (i,7 g]0,2], then there exists a constant c > such 
that, 



for Ixl < 1 . 



Proof. Since 



+00 



G}{x)< Pl{x)dt, 



the result follows from the well-known formula for Riesz 7-potential, see e.g. 
Landkof |La75j . □ 

Conclusion There exists a constant c > such that, if 7 < c?, 7 < 2, then: 



and 



dy I - |3,|<i-7 



if Ixl < 1 



if b|>l. 



Proof of Theorem 14. SI Now. we pass to the proof of the theorem and restrict 
to the case of b{t) = I, t > 0, that is, to stochastic heat equation. (Proof for 
the next two cases may be obtained in a similar way.) We have s{t; 7) 
t > 0, J > and therefore 

s{t;v{X)) ^e-"^^^', AeR'',i>0. 

By Theorem 14.71 if for some £ > and alH > 0, 



(ln(l + |A|)) 



fi{dX) < +00 , 



(4.17) 



then for all t > 0, the solution of the stochastic equation (in the general 
form), has a continuous version. Taking into account that a is a non-negative 
continuous function, one can replace (j4.17p by 
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(ln(l + |A|))i+-— i— Mc?A) (4.18) 
! l + i;(A) 

Since we have assumed that A = — (— zi)"/^, then v{X) = |A|", with a e]0, 2] 
and therefore (|4.18p becomes 

(ln(l + |A|))i+^— i— MdA)<+(X3. (4.19) 
1 + |A|" 

However, the condition (|4.19p holds for some e > if for some S > 

/i(c?A) < +0O . 



1 + |A| 



a — S 



In the same way as in the paper [KZOOj by Karczewska and Zabczyk, for some 
constant c > 0: 

where j := a — S. Taking into account Lemma 14.101 and Lemma 14.111 the 
result follows. □ 



4.2.4 Some special cases 

In this subsection we illustrate the main results obtained considering several 
special cases. 

Let us recall that the linear stochastic Volterra equation (|4.7p considered 
in the chapter has the following form 

X{t) = Xo+ f b{t-T)AX(T)dT + Wr{t), 
Jo 

where Xq e S'{M.'^), A is an operator given in the Fourier transform form 

T{AC)iX)^-v{X)TiOiX), Ce^'(K'), 

11 is a locally integrable function and Wr is an S"(R'^) -valued space homo- 
geneous Wiener process. This equation is determined by three objects: the 
spatial correlation F of the process Wr, the operator A and the function v 
or, equivalently, by the spectral measure /i, the function a and the function s, 
respectively. 

We apply our Theorems 14. 71 and 14.81 to several special cases corresponding 
to particular choices of functions v, a and of the measure /i. We will assume, 
for instance, that bit) ^ 1 or b{t) = t or b{t) = e"*, t > 0, that v{X) = |A|", 
a £]0, 2], A G R'^ and that the measure /i is cither finite or fi{dX) — j^dX, 
7 e]0,d[. Note that ii v{X) = |Ap, then A = A and if w(A) = |A|", a e]0,2[. 
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then A — — (— is the fractional Laplacian. In aU considered cases we 
assume that Hypothesis (H), on the function v, holds. 

Case 1 If (H) holds, the function a is given by (|4.5p and (|4.6p and the measure 
H is finite then TZ * Wr is a function- valued process. To see this note that by 
(H) and Theorem 14. 6[ the measure given by l|4.14p is finite. So, the result 
follows from Theorem |47 



Case 2 If (H) holds, the function a is given by (|4.5p and (|4.6p and fi is a 
measure such that for some e > 0, 

iln{l + \\\)^+' fi{d\) < +00, 

then for arbitrary t > 0, TZ* Wr{t) is a continuous random field. This follows 
immediately from Theorem 14.71 



Case 3 Assume that b{t) = 1 or b{t) = t or b{t) = e~\t>0,A = A (Laplace 
operator) and r{x) = rp{x) = (3 G [0,d[. Then function s is given by 

formulas (|4.10p . (|4.1ip and (|4.12p . respectively. Function v{\) = |Ap, and 
the spectral measure /x^ corresponding to Fp is of the form ^^{dX) = p^p^, 
with c/3 a positive constant. To simplify notation we assume that d > 2. 
Then TZ * Wr is a function- valued process if an only if /3 e]0, 2[, see ( [KZOOj '). 
Moreover, if /3 €]0, 2[ then for each t > 0, TZ * Wr{t), is a continuous random 
field. To prove this we use Theorem 14.81 and show that for some 6 > 0, 

T-T^ rf,ix)dx < +00 (4.20) 

|x|<l l-'-l 



and 

1 



-^T rf}{x)dx < +00 . (4.21) 



i.i>i 

Condition (|4.2ip is always satisfied because ()4.2ip is equivalent to: P > S — 2. 
Condition (|4.20p may be replaced by the following one: 

■ dx = c I , „ , c , r, r'^~^dr = c / , ^ dr < +oo , 



x|<l l-'-l JQ ' Jo ' 

equivalent to /3 < 2 — (5, which holds for sufficiently small S > 0. 
Case 4 Assume that b{t) — 1 and the operator A is given by the formula 

nAO{\) = -v{X)J^iO, 

where 

w(A) = {QX, X) + (1 - cos(A, x))v{dx) 
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and is a symmetric measure such that 

A 1) I'idx) < +oo . 



Then the equation (|4.7p has a function-valued solution if and only if 

fi{dX) < +00 . 



l + v{X) 



Additionally, if Xq = and 

(ln(l + |A|)i+^)— i--Ai(dA)<+oo, 

i l + V(A) 

then equation (|4.7[) has continuous version for each t > 0. 

In this situation, s(cr, w(A)) — e~'^^^'^\ By Theorem 14.71 the condition for 

function-valued solution of the equation (|4.7|) becomes: 



\Jo ) JM'' Jo 



and it is equivalent to 



, — — — - fi{dX) < +00 . 
l + w(A) 



4.3 Limit measure to stochastic Volterra equations 

This section is a natural continuation of the previous one. Description of 
asymptotic properties of solutions to stochastic evolution equations in finite 
dimensional spaces and Hilbert spaces is well-known and has been collected 
in the monograph [DZ96] . This problem has been studied for generalized 
Langevin equations in conuclear spaces also by Bojdecki and Jakubowski 
|BJ99] . The question of existence of invariant and limit measures in the space 
of distributions seems to be particularly interesting. Especially for stochastic 
Volterra equations, because this class of equations is not well-investigated. 

In the section we give necessary and sufficient conditions for the existence 
of a limit measure and describe all limit measures to the equation (j4.ip . Our 
results are in a sense analogous to those formulated for the finite-dimensional 
and Hilbert space cases obtained for stochastic evolution equations, see |DZ96[ 
Chapter 6]. 

Let us recall the stochastic Volterra equation (|4.ip in the simpler form 
((4J)) . that is, 

X{t)^Xo+ [ b{t-T)AX{T)dT + Wr{t) . 
Jo 
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As previously, we study this equation in the space 5"(R''), where e 5"(R'^), 
A is an operator given in the Fourier transform form (14. 8p . i.e., 

where is a locally integrable function and Wp is an 5" (R*^) -valued space 
homogeneous Wiener process. 

4.3.1 The main results 

In this subsection we formulate results providing the existence of a limit mea- 
sure and the form of any limit measure for the stochastic Volterra equation 
(|4.7p with the operator A given by (|4.8p . In our considerations we assume that 
the Hypothesis (H) holds. 

Let us recall the definition of weak convergence of probability measures 
defined on the space 5"(R'^) of tempered distributions. 



Definition 4.12. We say that a sequence {74}, t > 0, of probability measures 
on S"(R'^) converges weakly to probability measure 7 on 5"(R'') if for any 
function f G Cb(S") 



lim 

t — >+oc 



f{x)jt{dx) 



f{x)-f{dx) 



(4.22) 



More general definition on weak convergence of probability measures defined 
on topological spaces may be found, e.g. in |Bi68| or |KX95| . 

By we denote the law C{Z{t)) = A/'(0, Ft) of the process 



Let us define 



Z{t) 



Hoo{dX) := 



n{t- a)dWr{cr), t>0. 



(s(ct, a(A))) da 



fi{dX). 



(4.23) 



(4.24) 



Convergence of measures in the distribution sense is a special kind of weak 
convergence of measures. This means that 



ip{X)dfit{X) 



ip{X)dnooW 



(4.25) 



for any test function (p E 5(R'^). 

Now, we can formulate the following results. 
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Lemma 4.13. Let /i^ and fiao be measures defined by ^.14^ and ^.24^ , re- 
spectively. If ^J-Qo is a slowly increasing measure, then the measures lit /ioo; 
as t —f +00, in the distribution sense. 



Proof. First of all, let us notice that, by Theorem 14.61 Ht, t > 0, are spec- 
tral measures of stationary generalized Gaussian random fields. Moreover, the 
measures /it, t > 0, are slowly increasing. Since the function s(t, a(A)), t > 0, 
A G K'' is bounded, then the integral gt(A) = /q (s(t, a(A)))^ dr, for t < +00, 
is bounded, as well. In the proof we shall use the specific form of the measures 
Mt, i > 0, defined by gH]). 

We assume that the measure /^oo is slowly increasing, that is, there exists 
fc > 0: 



{l + lXrr'dii^iX)^ / (l + lAp) 



{s{t,a{X))fdT 



dn{X) < 



Hence, the function goc{X) = /g°°(s(T, a(A)))^ dr < 4«3 for /i - almost every A. 



In our case, because of formulae (|4.14p and (|4.24p . we have to prove the 
following convergence: 



lim 

t — ^+00 



^iX)gtiX)dfiiX) 



<y9(A).goo(A)d/i(A) , 



(4.26) 



where (p € S'(R''), and gt and ^oo are as above. 

In other words, the convergence (|4.25p of the measures Ht, t > 0, to the 
measure /ioo in the distribution sense, in our case is equivalent to the weak 
convergence (|4.26|) of functions gt, t > 0, to the function goo- 

Let us recall that the function s determining the measures fit,t > 0, and 
/Xcx3, satisfies the Volterra equation (|4.9p (see Hypothesis (H)): 



+ 7 / b{t~ T)s{T)dT = 1 



Additionally, by Lemma 2.1 from |CD97j . limt^+00 s(<) ~ 0. 
For any ip e S'(R'^) we have the following estimations 



^(A)gt(A)d/x(A) 



viX) gooiX)dniX) 



< 



< 



\^{X)\\gt{X)-goo{X)\dpi{X) = 



< 



l^(A)| 
l^(A)| 



t /• + 00 

2 / t^t^ „M^^^2 



{s{T,a{X))y dr ■ 



is{T,a{X))ydT 



d/i(A) < 



{s{T,a{X))y dr dfi{X) . 



(4.27) 
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The right hand side of (|4.27p tends to zero because 

/+CX; 
(s(T,a(A)))2dT 

tends to zero, as t +00. 

Hence, we have proved the convergence (|4.26p which is equivalent to the 
convergence (j4.25p of the measures /if, as t — > +00, to the measure /ioo in the 
distribution sense. □ 



Lemma 4.14. Let Ft, F^o be covariance kernels of the stochastic convolution 



I14-23 ) for t < +00 and t — +00, respectively, and let lit, fio^ be defined by 
14-14 ^ fl^*^ \4-'M^ - Assume that /ioo is a slowly increasing measure on . 
Then Ft — > i^oo, as t — )■ +00, in the distribution sense if and only if the 
measures fit ^ fJ-cxj, for t —>■ +00, in the distribution sense. 

Proof. The sufficiency comes from the convergence of measures in the dis- 
tribution sense which, in fact, is a type of weak convergence of measures. 
Actuahy, the convergence of /it, t > 0, to the measure /ioo in the distribu- 
tion sense means that (/if,(p) (/ioojip) for any ip G S'(M''). Particu- 
larly, because the Fourier transform acts from S{M.'^) into 5(R'^), we have 

{lJ,t,J^{(p)) (^fjL^,T{ip)) for any Lp £ 5(8*^). This is equivalent to the 

convergence {J^~^{fit),(p) *^±?° {T^^ {fioo) , v)) , ^ e 5(R'') . 

This means the convergence of the Fourier inverse transforms of considered 
measures /it, as t — *■ -l-oo, to the inverse transform of the measure /ioo in the 
distribution sense. 

Because the measures /it, t > 0, and fioo are positive, symmetric and slowly 
increasing on W^, then their Fourier inverse transforms define, by Bochner- 
Schwartz theorem, covariance kernels Ft — JF^^(/if), i > 0, and r'oo = 
^~^{tJ'oo)i respectively. Hence, Ft — s- -Too b& t +00, in the distribution 
sense. 

The necessity is the version of Levy-Cramer's theorem generalized for a 
sequence of slowly increasing measures {/it}, i > 0, and their Fourier inverse 
transforms which are their characteristic functionals. □ 

Now, we are able to formulate the main results of the section. 

Theorem 4.15. There exists the limit measure Vao = A/'(0, r'oo); ^he weak 
limit of the measures Vt = A/'(0, /^t); as t +00, if and only if the measure 
/ioo defined by \4-'M^ is slowly increasing. 

Theorem 4.16. Assume that the measure /ioo defined by ^.24^ is slowly in- 
creasing. Then any limit measure of the stochastic Volterra equation ^.7^ is 
of the form 

moo*AA(0,roo) , (4.28) 
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where rrioo is the limit measure for the deterministic version of the equation 
1^. 7| j with condition ^.8^ , and J\f{0, F^) is the limit measure of the measures 
Vf, as t +00 . 

We would like to emphasize that Theorems 14.151 and 14.161 have been for- 
mulated in the spirit analogous to well-known theorems giving invariant mea- 
sures for linear evolution equations, see e.g. [DZ96| or |BJ99| . Such results first 
give conditions for the existence of invariant measure and next describe all 
invariant measures provided they exist. Our theorems extend, in some sense. 
Theorem 6.2.1 from |DZ96| . Because we consider stochastic Volterra equations 
we can not study invariant measures but limit measures. 



4.3.2 Proofs of theorems 

Proof of Theorem 14.151 

{=>) Let us notice that, by Theorem 14.61 the laws vt = A/'(0, ri), t > 0, 
are laws of Gaussian, stationary, generalized random fields with the spectral 
measures /it and the covariances Ft- The weak convergence ()4.22|) is equiva- 
lent to the convergence of the characteristic functionals corresponding to the 
measures ft, i > and z^oo, respectively. Particularly 



hif) *^±?° Voo{^) for any Lp e 

We may use the specific form of the characteristic functionals of Gaussian 
fields. Namely, we have 

i>tM=Ee*<^W^^> =cxp(^-igt(^,^)) = exp (^-i(rt, ^ * 'Pc,))) , 

where t > 0, ip € S'(IR'*) and Z(t) is the stochastic convolution given by (|4.23p . 
Analogously 



i>^(<p) =cxp(^-^{Foo,'P*(p(s))^ , (p e 



Hence, we have the following convergence 

exp (^-^{Ft,(p* ipi^s))^ exp (^-^{Foc^p * (p(s)) 

for any ip e S{R'^). 

Because Ft, t >0, are positive-definite generalized functions then F^o is a 
positive-definite generalized function, too. So, by Bochner-Schwartz theorem, 
there exists a slowly increasing measure /^oo such that F^o — ^"^(Moo)- 

{^) Assume that the measure /Xqo, defined by the formula (|4.24p is 
slowly increasing. Then, by Bochner-Schwartz theorem, there exists a positive- 
definite distribution F^c on S such that -Too = ^~^(.l^oo) and 
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Now, we have to show, that Foe is the hmit, in the distribution sense, of 
the functional r^, < > 0. In order to do this, by Lemma [4.141 we have to 
prove the convergence of the spectral measures /if — *■ iioc, SiS t +oo, in 
the distribution sense. But, by Lemma [4. 131 the measures fit, t > 0, defined 
by l|4.14p . converge to the measure /Zqo in the distribution sense. This fact 
imphes, by Lemma [4.1 41 that Ft —^ F^c, as t — > +oo, in the distribution sense. 

Then, the following convergence 



holds for any (p e S'(R''). This means the convergence of characteristic func- 
tionals of the measures = -^{O-.Ft), t > 0, to the characteristic functional 
of the measure Voo = A/'(0, -Too)- Hence, there exists the weak limit v^o of the 
sequence vt, t > 0, and Voo — A/'(0, i^oo)- □ 

Proof of Theorem 14.161 

Consider a limit measure for the stochastic Volterra equation (|4.7p with the 
condition (|4.8p . This means that we study a limit distribution of the solution 
given by (|4.13p to the considered equation (|4.7p . 

Let us introduce the following notation for distributions, when < t < oo: 
r]t = C{X{t)) means the distribution of the solution X{t); nit — C(TZ{t)XQ) 
denotes the distribution of the part Tl{t)Xo of the solution X{t) and 
Pt = C{Z{t)) = TZ{t — T)dWr{T)) is, as earlier, the distribution of the 
stochastic convolution Z{t), that is, vt — N'{0, Ft). 

We assume that rjoa is any limit measure of the stochastic Volterra equation 
(|4.7p with the condition (|4.8p . This means that distributions rjt of the solution 
X{t), as t ^ +00, converge weakly to ijoo- 

We have to show the formula (|4.28p . that is the distribution 7700 has the 
form rjoo = Woo * A/'(0, Too ) . 

The distribution of the solution (|4.13p can be written 



for any <t < +00. 

Because the initial value Xq is independent of the process Wr (t) , we have 





C{X{t))^£{TZ{t)Xo)*£{Z{t)) 



or, using the above notation 



r]f = rrit * 1/4, for any < t < +00 . 
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This formula can be rewritten in terms of characteristic functionals of the 
above distributions: 

fit{ip) = rhtiip) i>ti^) , (4.29) 

where ip e S{M.'^) and <t < +00. 

Then, letting in (|4.29p t to tend to +00, we have 

where f]oo{^) is the characteristic functional of the limit distribution 770^, 
C{(f) = limj^oo TOt((^) and Voai^) is the characteristic functional of the limit 
measure Voo = 7V(0, Too); moreover Voo{<p) = cxp(-i(roo, * y^is)))- 

Now, we have to prove that C{ip) is the characteristic functional of the 
weak limit measure moo of the distributions mt = C{TZ{t)Xo)- 

In fact, 

C(¥') =??oo((/')exp Q(r'oo,'/3*¥'(s))^ , 

where the right hand side of this formula, as the product of characteristic 
functionals, satisfies conditions of the generalized Bochner's theorem (see e.g. 
|It84| ). So, using the generalized Bochner's theorem once again, there exists 
a measure rrioo in S'{M.'^), such that C{(p) — rhao, as required. Hence, we have 
obtained rjao = * A/'(0, Foe)- □ 



4.3.3 Some special case 

Stochastic Volterra equations have been considered by several authors, see 
e.g. |CD96l I^DQTl ICD00j . [CDP97] and [KZOOaj . and are st udied in connec- 
tion with problems arising in viscoelasticity. Particularly, in |CD97| the heat 
equation in materials with memory is treated. In that paper the authors con- 
sider an auxiliary equation of the form 

z{t)+ f [iuc{t-T)+l3{t~T)]z(T)dT = l, (4.30) 
Jo 

t > 0, where /i is a positive constant and c, P are some functions specified 
below. 

Let us notice that if in the Volterra equation (14. 9p we take 6(r) = ^ [/i c(r)-f 
P{t)], we arrive at the equation (|4.30p . On the contrary, if we assume in the 
equation (|4.30p that /3(t) = 0, /x = 7 and &(t) = c(t), we obtain the equation 

(|tgD . 

Assume, as in |CD97] . the following HYPOTHESIS (HI): 

1. Function /3 is nonnegative nonincreasing and integrable on R"*". 

2. The constants /i, co are positive. 
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3. There exists a function S G L^(M+) such that: 

y+OO 

c{t) := Co 



/-t-oo 

/ S{a)d(T and Cqo := cq — / |(5(CT)|(iCT > 
Jo Jo 



Proposition 4.17. (jCDQTf, Lemma 2.1) Let functions f3, 6 and c be as in 
Hypothesis (HI). Then the solution to ^J^SU^ satisfies: 

1.0<\z{t)\<l, t>0; 

2. /+°°|z(t)|dt< (^coo)-i <+oo. 

In the next result we will use the above assumption and Proposition [HIT] 
of Clement and Da Prato and follow the spirit of their argumantation. 

Proposition 4.18. Assume that the stochastic Volterra equation 7| ) has 
the kernel function v given in the form 

b{t)^co-[ \S{a)\da>0, cq > where S E L^{R+) , 
Jo 

and the operator A is given by ^4-^ - ihis case the limit measure given 
by the formula a slowly increasing measure. 

Proof. This Proposition is the direct consequence of the definition (|4.24p of 
the measure /ioo and Proposition l4.17l In fact, from Proposition 14. 1 71 we have 

|s(T,7)|dT< (7Coo)-\ (4.31) 







where 7 satisfies Hypothesis (HI), so Coo is finite. Hence, the right hand side of 
(|4.3ip is finite for any finite 7. In our case, because A satisfies (|4.8p . 7 — a{\). 
From the definition ()4.24p of the measure /ioo we have: 



(l + |An-'=d/.oo(A)= / (l + |Ap) 



{s{t,a{\))f dr 



dn{X) (4.32) 



for k > 0. Let us notice that, by Proposition I4.17| < |s(t, a(A))| < 1 for 
t > 0. So, (s(i,a(A)))2 < |s(t,a(A))|. Therefore, because (|i3T|) holds and the 
measure ^ is slowly increasing, the right hand side of (|4.32p is finite. Hence, 
the measure /ioo is slowly increasing, too. □ 



4.4 Regularity of solutions to equations with infinite 
delay 

4.4.1 Introduction and setting the problem 

In this section we consider the following integro-differential stochastic equation 
with infinite delay 
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X{t,0)^ f b{t- s)[AX{s,e) + Wr{s,9)]ds, t > 0, deT'^, (4.33) 

where b G L^(R+), A is the Laplace operator and T'' is the d-dimensional 
torus. In (I4.33p . Wr is a spatially homogeneous Wiener process with the space 
covariance F taking values in the space of tempered distributions S'{T'^) and 
Wr denotes its partial derivative with respect to the first argument in the 
sense of distributions. Such equation arises, in the deterministic case, in the 
study of heat flow in materials of fading memory type (see |CD88] . [Nu71| ). 

In this section we address the following question: under what conditions 
on the covariance F the process X takes values in a Sobolev space H°'{F'^), 
particularly in L'^{F'^)7 

We remark that the knowledge of the regularity of solutions is important 
in the study of nonlinear stochastic equations (see e.g. |DF98] and |MS99| 1. 

We study a particular case of weak solutions under the basis of an explicit 
representation of the solution to (I4.33P (cf. Definition I4.20|) . 

Observe that equation (I4.33|) can be viewed as the limiting equation for 
the stochastic Volterra equation 



X{t,9)= f b{t-s)[AX{s,9) + Wris,0)]ds, t>0, 9eF'^. (4.34) 

"'0 

If b is sufficiently regular, we get, by differentiating (|4.34p with respect to t, 



dX I 

— it,e)^biO)[AX{t,9) + Writ,9)]+ J b'{t - s)[AX{s,9) + Wr{s,9)]ds, 

(4.35) 

where i > and 6* € T''. 

Taking in ((i35)l . b{t) = 1 we obtain 
dX 

— {t,9) = AX{t,9) + Wr{t,9), t>0,9eT^, ^^^^^ 
X{0,9)^0, 9(zF'^. 

Similarly, taking b{t) = t and differentiating ()4.34p twice with respect to t we 
obtain 

^{t,9) = AX{t,9) + Writ,9), t>O,0GT^ 

X(0,6i) = 0, 61 gT-^, 

^(0,6») = 0, 9eT''. 

It has been shown in [KZOO| Theorem 5.1] (see also }KZ01|, Theorem 1]) that 
equations and (|07)) on the d-dimensional torus T'^ have an H°'+^{F'^)- 

valued solutions if and only if the Fourier coefficients (7„) of the space covari- 
ance F of the process Wp satisfy 
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7n(i + H'r <00. 



(4.38) 



Observe that for both, stochastic heat (|4.36p and wave (|4.37p equations, 
the conditions are exactly the same, despite of the different nature of the 
equations. On the other hand, the obtained characterization form a natural 
framework in which nonlinear heat and wave equations can be studied. 

In this section, we will prove that condition (I4.38P even characterizes 
-ff"+^(r'')-valued solutions for the stochastic Volterra equation l|4.33p . pro- 
vided certain conditions on the kernel b are satisfied. This is a strong contrast 
with the deterministic case, where regularity of (|4.33|) is dependent on the 
kernel b. The conditions that we impose on b are satisfied by a large class of 
functions. Moreover, the important example b(t) = e^* is shown to satisfy our 
assumptions. 

We use, instead of resolvent families, a direct approach to the equation 
(|4.33p finding an explicit expression for the solution in terms of the kernel 
6. This approach reduces the considered problem to questions in harmonic 
analysis and lead us with a complete answer. 

Let {f2,!F,{!Ft)t>o,P) be a complete filtered probability space. By T"* 
we denote the c?-dimensional torus which can be identified with the prod- 
uct (— 7r,7r]'*. Let D{T'^) and D'{T'^) denote, respectively, the space of test 
functions on T"^ and the space of distributions. By (^, (j)) we denote the value 
of a distribution ^ on a test function. We assume that Wp is a Z?' (r'')-valued 
spatially homogeneous Wiener process with covariance F which is a positive- 
definite distribution. 

As we have already written, any arbitrary spatially homogeneous Wiener 
process Wr is uniquely determined by its covariance F according to the for- 
mula 



homogeneous process, the distribution F — F{6 — rj) for 6, -q G F'^. 

The space covariance F , like distribution in D'{F'^), may be uniquely 
expanded (see e.g. |GW99| or [Sc65| ) into its Fourier series (with parameter 
w = 1 because the period is 27r) 



convergent in D'{F'^). In (|440l) . (n, 0) = Y,t=i niOi and Z'^ denotes the product 
of integers. 

The coefficients 7„, in the Fourier series (|4.40p . satisfy: 

1. 7„ = 7_„ for n eZ'^, 

2. the sequence (7„) is slowly increasing, that is 






(4.40) 




for some r > 0. 



(4.41) 
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Let us introduce, by induction, the following set of indexes. Denote 
Zl := N, the set of natural numbers and define Zf+i := (Z^ x Z'') U {(0, n) : 
n e Zf}. Let us notice that Z'^ = Zf U (-Zf) U {0}. For instance, for d = 2, 
Z2 = NxZU{(0,n) :nGZ}. 

Now, the spatially homogeneous Wiener process Wr corresponding to the 
covariance F given by (|4.40p may be represented in the form 



Wr{t,e) = Vl^f3o{t) + y2^[cos(n,0)/3i(i) + sin(n,0)/32(t)] ,(4.42) 



In (|4.42|) Pq, Pli, P'^,n d Zf , are independent real Brownian motions and 70, 79 



are coefficients of the series (|4.40p . The series (j4.42p is convergent in the sense 
of D'{T'^). 



Because any periodic distribution with positive period is a tempered dis- 
tribution (see, e.g [GW99j ). we may restrict our considerations to the space 
5"(r'*) of tempered distributions. By S'(T'') we denote the space of infinitely 
differentiable rapidly decreasing functions on the torus T''. 

Let us denote by i?" = H°'{T'^), a e M, the real Sobolev space of order a 
on the torus T'^. The norms in such spaces may be expressed in terms of the 
Fourier coefficients (see |Ad75j ) 



where ^„ = + i^i ^=C-n, n£'L<^. 

There is another possibility to define the Sobolev spaces (see, e.g. |RS75| ). 
We say that a distribution ^ g S'{T'^) belongs to H"',a G K, if its Fourier 
transform ^ is a measurable function and 



4.4.2 Main results 

If 6 G L[q^(K+) and /x g C, we shall denote by r(t, fi) the unique solution in 
L\^^{W+) to the linear Volterra equation 



t > and 61 e T'^. 





r{t, fi) b{t) + n [ b{t- s)r{s,n)ds, t>0. (4.43) 
Jo 
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In many cases the function r(t, /i) may be found expHcitly. For instance: 



b{t) =1, 




b{t) ^t, 




b{t) ^ e-\ 




b{t) = te~*, 





sinh t 

e(i+M)* 



_ smh^t 
e ^— , MT^O. 

For more examples, see monograph |Pr93| by Priiss. 

Let us denote by f{k),k G Z, the fc-th Fourier coefficient of an integrable 
function /: 



Jo 

Given b e L'^(R+), we find that, for F{t) J*^^ b{t - s)f{s)ds, we have 

F{k) = b{ik)f{k), keZ, (4.44) 

/"OO 

where 6(A) := / e"'^* b{t)dt denotes the Laplace transform of b. 

In what follows we will assume that b{ik) exists for all fc S Z and suppose 
that A 6(A) admits an analytical extension to a sector containing the imag- 
inary axis, and still denote this extension by 6. We introduce the following 
definition. 

Definition 4.19. We say that a kernel 6 G L^(]R+) is admissible for equa- 
tion (im if 

lim \nf / [r{s,-\nf )fds =: Cb 

exists. 
Examples 

1. In the case b{t) = e~* we obtain 



and hence Cf, = |. 
2. In the case b(t) ~ te~* we obtain 

and hence Ch — \- 
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Denote by Wn{t, 9) cos(n, 0) f3l{t) + sin(n, 0) l3l{t), n e , that is the 
n-th element in the expansion (|4.42p . 

Definition 4.20. By a solution X(t,0) to the stochastic Volterra equation 
i4.33\ ) we will understand the process of the form 



X{t,0)^^Po{t)+Y. f Ki-s, 



where the function r is as above, t > and G T*^. 



\n\^)dWnis,0), (4.45) 



The process X given by (|4.45p is a particular form of the weak solution to 
the equation (|4.33p (cf. [KZOOj ) and takes values in the space S"(T^). 

The following is our main result. 

Theorem 4.21. Assumeb E L^(R+) is admissible for 1^4. Then, the equa- 
tion has an H°'^^{T'^)-valued solution if and only if the Fourier coeffi- 
cients (7„) of the covariance F satisfy 

^ 7„(1 + Inp)" < oo. (4.46) 



In particular, equation has an F'^{T'^) -valued solution if and only if 



In 

: — TTT < oo. 



Proof. We shall use the representation (|4.42p for the Wiener process Wr{t, 0) 
with respect to the basis: 1, cos(n, 6*), sin(n, 0), where n g Zf and E T'^. 
Equation (|4.33p may be solved coordinatewise as follows. 
Assume that 

X{t, 0)^Y^ [cos(n, 61) Xl{t) + sin(7i, 0) Xl{t)] + Xo{t). (4.47) 

Introducing (|07l) into (|433l) . we obtain 
cos(7i, 0)Xl (t) + sm{n, 0)Xl (t) 

= -|n|2 / b{t-s)[coa{n,0)X^{s)+sinin,0)X^{s)]ds 



+ v/2^ / h{t - s)[cos{n,9) Pl^is) + sm{n,0) f3lis)]ds, 

or, equivalently 
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[cos(n, 9), sin(n, 



'Kit) 

Kit) 



= -|np / s)[cos(n,6'),sin(n,0)] 



+ \/ 27n / 6(i — s)[cos(n, 0), sin(n, 0)] 



ds 



Denoting 

we arrive at the equation 



and f3n{t) 



X„it) = -\n\^ b{t - s)Xn{s)ds + b{t ~ s)dpn{s) . (4.48) 

J — oo J —oo 

Taking Fourier transform in i, and making use of (|4.43p with ji = — |np and 
(|4.44p . we get the following solution to the equation (|4.48p : 



r{t- s,-\n\^)^/2^d(3n{s) ^ / r{s, -\n\^)y^d/3n{t - s) 



Xn{t) 



Hence, we deduce the following explicit formula for the solution to the equation 
(14331) : 



cos{n,0) / r(s, -|n|^)d/3^(t - s) 



+ smin,9)l ris,~\n\')dp'^{t^s)^ 

(4.49) 

Since the series defining the process X converges in S"(T'^), P-almost surely, 
it follows from the definition of the space iJ" that X{t) E P-almost 
surely if and only if 



^ (1 + InH^+^T. 



t X 2 

r{t-s,-\n\')dpl,{s) 



r{t-s,-\n\')dPl{s) 



(4.50) 



< oo. 



Because the stochastic integrals in (j4.50p are independent Gaussian random 
variables, we obtain that (j4.50p holds P-almost surely if and only if 



J2a + \n\Y^'ln E 



X 2 

rit-s,-\n\^)d/3lis)] (4.51) 



t X 2- 

+ ( / r{t-s,-\n\')dPl{s) 



< oo. 
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Or equivalently, using properties of stochastic integrals, if and only if 

^ (1 + |nn"+i7„ / [r{s-\n\'')fds<^. (4.52) 

Since h is admissible for the equation (|4.33p . we conclude that (I4.52p holds if 
and only if 



and the proof is completed. □ 

Concerning uniqueness, we have the following result. 

Proposition 4.22. Assume b e L^(R+) is admissible for Ii4.33\ ) and the fol- 
lowing conditions hold: 

(l) ^ 7n(l + IH'r<00, 

(ii) C C \ : n e 'L''}. 

Then, {4 -33^1 has a unique H°'^^{T'^)-valued solution. 
Proof. Let X{t,6) be solution of 

X{t,e)^ f b{t- s)AX{s,9)ds. 

J — oo 

Taking Fourier transform in 9 and denoting by Xn{t) the n-th Fourier coeffi- 
cient of X{t, 9){t fixed), we obtain 

Xn{t)^~\n\^ f b{t - s)Xn{s)ds 



for aU n e Z'^. Taking now Fourier transform in t, we get that the Fourier 
coefficients of Xn (t) ( n fixed) satisfy 

(1 + \n\%{ik))Xn{k) = 

for all fc e Z. According to (ii) we obtain that X„ {k) = for all fc e Z and all 
n G Z'*. Hence, the assertion follows by uniqueness of the Fourier transform. 
□ 

The following corollaries are an immediate consequence of Theorem 14.211 
The arguments are the same as in |KZ01] . We give here the proof for the sake 
of completeness. 

Corollary 4.23. Suppose b G L^(M_|_) is admissible for ^.33^ and assume 
r G L'^[T'^). Then the integro- differential stochastic equation ji4-33\ ) has a 
solution with values in L'^iT'^) for d = 1, 2, 3. 
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Proof. We have to check equation (|4.46p with a — —1. Note, that if e 
L2(T'') then r = (7„) e P(Z''). Consequently 

V < f V -7 V ^- 

^ l + |n|2 - ^ ^Z^ (1 + |„|2)2^ ■ 

But Z^nez'i < and 

7 1 inx < oo if and only if 2q > d. (4.53) 

Hence, the result follows. □ 

Corollary 4.24. Suppose b G L^(M-|_) is admissible for ^.33^ and assume 
r G ^{lA) for 1 < p < 2. Then the integro- differential stochastic equation 
{4-33^ has a solution with values in L'^{T'^) for all d < ^zj- 

Proof. Note that 

1+|„|2 - ^2^ T» ^ ^2^ (H_|„|2)gi ' 

where ^ + ^ = 1. Hence the result follows from (|4.53|) with q — □ 

For a = — 1, the condition (j4.46p can be written as follows. 

Theorem 4.25. Let b be admissible for the equation {4-33^ - Assume that the 
covariance F is not only a positive definite distribution but is also a non- 
negative measure. Then the equation Ii4.33\ ) has L^{T'^) -valued solution if and 
only if 

(r,Gd)<+(X), (4.54) 

where 

Gd(a;)= V / e-*e- « dt, x£T<^. (4.55) 

i^.io y(W 

The proof of Theorem 14.251 is the same that for Theorem 2, part 2) in 
|KZ01| . so we omit it. For more details concerning the function Gd we refer 
to [KZOT] and prf5] . 

Additionally, from properties of function Gd defined by (|4.55p and the con- 
dition (|4.54p we obtain the following result (see |KZOO[ Theorem 6.1]). 

Corollary 4.26. Assume that F is a non-negative measure and b is admissi- 
ble. Then equation \4-3!^ has function valued solutions: 

i) for all F if d ^ 1; 

ii) for exactly those F for which J|g|<]^ In \9\ F{d9) < +oo if d = 2; 
Hi) for exactly those F for which J^g^^i \e\^-^ F{dO) < +oo if d> 3. 
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In what follows, we will see that formula (|4.49p also provides holderianity 
of X with respect to t. In order to do that, we need assumptions very similar 
to those in [UD96j . 

HYPOTHESIS (H2) Assume that there exist 5 € (0, 1) and Cs > such that, 
for all s G {—oo,t) we have: 



(^) j [r{t ~ T,~\nf)fdT < Cs\nf^'-'^ \t ~ sf ; 

(n) r [r{t-T,-\n\^)-r{s-T,^\n\^)]^dT <Cs\n\^'^^-^^\t- s\^. 

J —oo 

Proposition 4.27. Assume that ^ ^" < +oo. Under Hypothesis (H2), 

the trajectories of the solution X to the equation ^.33^ are almost surely rj- 
Holder continuous with respect to t, for every rj G {0,6/2). 

Proof. From the expansion (|4.49p and properties of stochastic integral, we 
have 

E \\X{t,9) - X{s,0)\\l, ^E\\^{(3f,{t) - Po{s)) 



cos{n,9)[ r{t-T-\n\')d(3l,{T)- r{s-T-\n\')d(il,{T) 



+ sin(n,0) / r{t-T-\n\')dPi{T)- / r{s-T-\n\')dPi(T) 

\j—oo J —OO 

{27ry{lo\t-s\ 

f ^ 27„ f [r(t-r,Hn|2) - r{s-T,-\nf)fdT+ f r\t~T,-\n\^)d7 

^riT^ lJ—OO Js 



|2 

Il2 



According to assumptions (i) and (ii) of the Hypothesis (H2), we have 
E\\X{t,d)-X{s,9)\\l.<Cs ^ 27„|np(*-i)|i-s|^ 

Because X is a Gaussian process, then for any m Cz N , there exists a constant 
Cm > that 



E||X(t,0)-X(s,0)||i? <C„ 



t~S 



Taking m such that m6 > 1 and using the Kolmogorov test, we see that the 
solution X{t, 0) is T^-Holder continuous, with respect to t, for rj = 5/2— l/(2m). 
□ 



92 4 Stochastic Volterra equations in spaces of distributions 

Example Let us consider the particular case b{t) = e~*, t > 0. Then, by 
previous considerations, r(t, — |np) = e*^^~l"l One can check that in this 
case the Hypothesis (H2) is fulfilled. 

Remark We observe that the condition (i) in Hypothesis (H2) is the same as 

|np [\r{s,-\n\')rds<Cs\n\^'\t\' 
Jo 

and hence it is nearly equivalent to say that the function b is admissible. 
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